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ON TWO-DIMENSIONAL • FLOWS Of COMPRESSIBLE FLUID-S 

By' Stefan Bergman 


SUMMARY 


This report is devoted to the study of two-dimensional 
steady motion of a compressible fluid. 


It is shown that the -complete flow pattern around a 
closed obstacle cannot be obtained by the method of Chaplygin. 
In order to overcome this difficulty, a formula for the 
stream- function of a two-dimensional subsonic flow is derived. 
The formula involves an arbitrary function of a complex vari- 
able and yields all possible subsonic flow patterns of certain 

r — i 


types. It is a generalization of the expression Im 
for the stream function of an incompressible fluid, 
is the velocity vector and g an arbitrary analytic 


[*<n] 

(Here V 
funct i on . 


) 


Conditions are given so that the flow pattern in the 
physical plane will represent a flow around a closed curve. 

The formula obtained can be employed for the approximate 
determination of a subeohic flow around an obstacle. The 
method can be extended to partially supersonic flows. 


INTRODUCTION 


The theory of irrotational t wo-dimen'si onal flows of an 
incompressible fluid is based on the theory of analytic func- 
tions of a complex variable. 

The relation between these two theories is given by the 
fact that the stream function \{f(x,y) of flow satisfies the 

Laplace equation :(3 s \jr/0x 2 ) + (S s ^/3y s ) = 0. Hence the imag-r 
inary part of an analytic function f(x + iy) is a stream 
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function of a possi'Dle flow, and all flow patterns can be ob- 
tained in this way. 

For certain purposes, however, it 1 b useful to modify 
this approach. The stream function may be considered as a 
function of the components and v 3 of the velocity vec- 
tor Again \J/ satisfies the Laplace equation 

(d 2 \(//dv 1 a ) + (6 s \Jf/8v g 3 ) = 0. Therefore, it is possible to 
choose as ^(v^Vg) the imaginary part of an analytic func- 
tion g(v), v being a complex variable in the (tj , v g )- 
plane. In this way the flow pattern in the ( v x , v £ ) -plane 
(hodograph. plane) is obtained. In order to find the real 
shape of the streamlines it is neceEBary to derive from 
Im g(v) the corresponding function., of x and y. This 
transition does not involve any serious theoretical diffi- 
cult i e s , 

In the case of a potential flow of a compressible fluid 
the first method (construction of the flow pattern directly 
in. the. physical plane) leads to a rather complicated nonlin- 

1 ear partial differential equation. The second approach (con- 
struction of the flow pattern in the hodograph plane) reduces to 
the ' integf ati on of a linear partial differential equation. 

'( See Chaplygin, reference 1.) Hence, the use of the hodo- 
graph method permits the application of various results from 
the, theory of linear part ial' di f f er ent ial equations." For 
instance, a stream function in the hodograph plane can be ob- 
tained as a linear combination of particular solutions of the 
linear equation mentioned. Chaplygin was the first to con- 
struct a set of siich particular solutions'. Two other methods 
of constructing such sets have -'been giveh by the present', 
author. (See reference- 2, pp . 16-20 and 23-24, and reference 
3 , sec . 2 . ) 

However, Chaplygin's method and both methods given in 
references 2 and 3 are not satisfactory in one respect. In 
general, the stream-function will be represented by an infi- 
nite series of particular solutions, and such a series will 
converge only within a part of the domain in which the flow 
is defined 1 , 


X A hodograph of a flow around a profile is (in general) 
a multiply covered domain (see fig. lb and 2b) the branch 
points of which are not necessarily located either at the 
or-i-gin or at infinity; on the other (continued on next page) 
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To obt-ain results pertaining to the actual flow, a, rep-, 
resentation of the stream function as a whole is indisp'ensi- 
ble . A representation fulfilling these requirements is given 
in this paper. (See also reference 2, sec. *6, and reference 
3, sec, 4.) 


If a linear relation "between the pressure p and. the 
specific volume l/p is assumed: 

p = A + cr/p (l) 

(A,cr constants), then the hodograph equation. coincides with 
the Laplace equation. Assuming relation (l) and using the 
theory of functions of a complex variable, Von Karman (refer- 
ence 4) and Tsien (reference 5) obtained the compressible 
flow past an elliptic cylinder. Equation (l) is a very tfough 
approximation to the actual pr e s sur e -den sity relation and can 
be used only in cases where the local velocity is far below 
that of sound. 

In the present report a general pressure-density relation 


p = A + ap k (2) 

is used (A, <j, fc are constants). (Equation ( 2 ) contains as 
a special case the adiabatic relation p = crp 1, .) Assuming 
(2) gives a general formula for' the stream function. This 
formula expresses the stream ..function of a compressible flow 
in terms of an arbitrary analytic function of a complex vari- 
able. 


The representation obtained is, in general, valid in the 
whole region where the flow is subsonic- and- in some cases can 
be extended into a supersonic region also. 

This investigation, conducted at the Brown University 
was sponsored by,, and conducted with the financial assistance 
of the National Advisory Committee for Aeronautics. 

(continued from page 2) hand, the Chaplygin solutions 
yield flows which (in the hodograph plane) either are single- 
valued, or multi-valued with a branch point at the origin or 
at infinity. In order to represent such flow patterns, sev- 
eral series development, each of which represents the. stream 
function \J/ under consideration in a certain part of the 
domain In which is defined, is needed. 
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I take the opportunity to express my gratitude to Mr. 
Leonard Greenstone for his assistance in the preparation of 
the present paper. 


NOTATION- 


Remark : In dealing with differential equations, the following 

complex notation is often used: 


_ Bu - f 
8.Z 2 \ 3x 


i 


au\ 

dy/' 


z dz 2 \ dx 9y/ 


U.J = I a if _Sfu. ' + 

ZZ 4 4Vs X 2 3y®/ 


z = x + iy, z = x - i y 


a = 


a o ~ ~ 1 ' ) 


’ • ■'! . 

v 3 p 


speed of sound; (equation (28)) 


a 0 speed of sound at a stagnation point 

c (See ( 94) ) 
c^ n ) (See (94).) 


exn(x) = e x ; e base of Naperian logarithms 


f (z) an arbitrary analytic function df • the complex vari- 
able z ‘ ; 

fV ■=• 3f/3t ; f t = df /at 

; • • ■ ’ • ; 
g constant of gravity 

,g(|j)' : 'an analytic function of the complex var iab.le . £; the 
■'result of applying ' ,i he transf or ma tion' z = z (£ ) to 
: f f ‘ ’ f(z) : ' >;r 

e {o) u) = scv ' 

g ^(0 (See ( 11-3 ) , f f .,); ; v . ■ 



HACA TU; Ko.r -972 


5 


h = ' f — for - k '>' .A o c c'a e i.o ha IT y • the boundary curve of 

J k + 1 a 'd&mdin. H. • 

k ratio of specific beat at' 'constant' pressure to constant 
volume 

l(H) =(^) 3 = U - M 2 (H))j (4.5). . 

\dES \p (H ) / 

/ «■ 

p pressure 

p Q pressure 11 at rest" 

r polar coordinate in the physical plane- 


s = 1 - J 1 - M 2 ; (151) 

* « 

v 

schlicht = univalent 

v sp’eed; magnitude of also,’ occasionally, the reduced 

speed v/a 0 

— ^ 

v x ,v 3 Cartesian components of V 

W st <J> + i\jf ’ 

x,y Cartesian coordinates in' the physical plane 

v r-, •y.cr- 

z = x + iy 
"z = x iy 

A constant in the pressure-density relation (22) (See 


also sec, 3,) 

D = = d± _ 3£_ - ’ (132) ■■ ■ 

3 ( x , y ) 3 x 3 y 3 y 3 x 


, v / 


E (See theorem (53).) Eg = Eg = §-5.- 


. . J + L ...\- 

E*. = exp! /kd ( g*I) j E; (69) 


3 H 


\-co 


30 ' ’ 

i 
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F = -(S. + N 3 ) = .fL.t.. L)..^ 4 .-[- (3k - 1 ) M 4 —4(3 - 2k)M 2 + 1 6~I ; ( 7 1 ) 

£ 64(1 - M 2 ) L J 

m 

F m (2A) = y e sXk ; ( l emma (67)) 

k=i _ 

F* (See theorem (85) . ) 

F a ( See (115), ff . ) 

G-(E) (See (53); also (124).) 

& 1 (E) (See (68).) 


ff s (E) (See (75).) 



(ill) 1 ; occasionally a domain in the 

holograph plane with boundary curve h 


H (See (115), f f . ) 


heat content 


i/g pressure he ad 


( top 11 - 1 \ 

\^g(k - 1 ) J 


Im the imaginary part of 

K kernel function (See reference 9; also sec, 3.) 
Kg ker.'nel' 'fund ti-o'rf of H;' r = P Kh(V,A)<1V 


L(\]/) = \|/pr + F\|/; (70) 

I 0 (\J/) + E(2?0 \|/^ + ^ (46) 

L m (\J/) = \J/^ + V; (74), ff, : 

L a . (See (115).) 


ft 
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Re the real part of 
Sfihlloht = lirtivdlsnt 
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V a ve = v, + iv_ velocity vector 
~f = ve- i6 = v x - iv 2 = *£ 


Vq, magnitude of the velocity at infinity 


W,W* (See (119).) 


2-e 


{’(Ktttt) ' 1s “ 


a, the angle which a doublet makes frith the real axis. (14) 

(Also, a real translation of the axis (See (ll5) ante.) 

e, aeB “a is a member of B" or "a belongs to B M 

£ = X ± i0 in which case ( = A ; ifi (This variable use of 
l merely means a reflection with respect to 
the real axi s . ) 

-S’ 

Q the angle which 7 makes with the real axis 


A (M) = A(v); (48), (49) 

p density; p = p^ £l -• 

p the density "at rest" 


tfc T ' 1) y 2 

2a® 


JET, 


(25) 


<J a constant in pressure-density relation p = A + CTp lc ; (22) 

potential function 1 ; also, the polar angle in the physical 
plane (polar coordinates) * “ 

stream function 1 

( p P 'p ' ' 

it* = expi / Hd^+^)U ; (69) I 

V / "•;■•;! " ' ' z"_ 

T circulation; in part II the Gamma fuhctiori r(x) = / e” 1: t x ‘' 1 dt ; 

(flee sec. 11, ff)< J 

o 


See r eaark-'i . 
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A_ Laplace operator: + — % = 4 ( 

dx® By® ^BzBi^ 

i • *, , , 

A(H) (See (45).) 

A s •= 3A/3H, Ag = 3A/30 

$ Potential Junction* (See r eaark 1 . ) 

Stream Junction (See r emar k 1 . ) 


Remark 1 : In the following, the potential and the stream 

functions <J>(x,jr) and xj/^x.y) (as well as several other 
variables whiqh ; ar« indicated in this section) are considered 
as functions of different pairs of variables. In passing 
from the physical to some other plane, new symbols should be 
introduced for 4> and \J/, since in different planes 4> and 
\j/ are different functions ,of their respective arguments, 
j or -instance , tpfa ? S‘sI ; hg' to the (v,0)-plane yields 


■ 4> ^ 1 ^ ' (v j 0 ) = d 
^ (i) ( v , 0 ) 


x(v,0), y(v,8) 


x ( v , 0 ) , y(v,e) 


.Jor the ■ sake -of brevity the author omits the superscript and 
always writes <$ and , no matter in which plane he con- 
siders these functions. 

• . POES WORD .... ■ 


The Itreaa function of an incompressible fluid flow 

is a solut i On ‘ of' "‘the Laplace equation 




' ^ - o •• 


Bx __ By 

There -efi'sts 'a general formula 




= I m 


f ( s) 


s = log v ' - "i0” 


■U) 


(4) 
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in terms of an arbitrary function f of one variable, for 
solutions of this equation. Here v is the speed and 0 
the angle which the velocity vector forms with the positive 
x-axis. 

In the case of a flow of a compressible fluid, the 
stream function is a solution of the system of equations 1 



where p Q , a Q , and k are constant. For air, k 


1.4. 


A generalization of formula (4) to the case of subsonic 
flows of a Compressible fluid is given in this paper.. Let 



be the local Mach number. 


If it 


is assumed that the flow is subsonic and that k = 1.4, func- 
tions A(M), H(M), Q,( n '(M), n = 1,2, . . . are determined 

(see table lb) so that for the solutions of (5) there Is ob- 
tained a representation 8 


1 St r cam fun ction and density p have t o be c on s i d- 

ered as iinknown in system (5). The elimination of p, In 
order to obtain one equation for is impossible. 


s For many purposes in evaluating (6), it suffices to 
take only a few terms of -the series. There also exist methods 
for improving the convergence of (6). As will be shown else- 
where, It is possible under rather general assumptions to 


interchange the 


lim 


m 


and the summation 


y 

n="i 


and thus 


obtain a new formula for \J/. In many instances, however, the 
formula in the original form is more suitable for applica- 


tions, since by a suitable choice of the 
to achieve faster convergence. 


m's it is possible 
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' \|f(Mve : >" = ’lim; Im (h(M) (£) 

' • ' ■ '■ m — 3> eo L L 


n=l ’ o o 




j (6) 


i = mm) - 


iG , x C m) s I log 
6 


1 - (1 - M 3 ) 1 + h(l - M a ) 




u 


1 + (l - M 3 )^ \ 1 - h(l 


- m 3 ) : 


-1^ '• --'V- - 


h = , k'> : i : 

^k + 1/ 


in terms of an arbitrary function of one variable 1 ... ’Since f 
the transition to the variables x, y does not involve any 
essential difficulty, (4) and (6). yield patterns for possible 
incompressible and subsonic (compressible) fluid flows. 

' * 

Formula (6) is of interest not merely as a tool for com- 
puting ’examples of flows cf -a ’compressible, fluid*, but -it- may- 
be 'consider ed also as an - operat i on ' which t rahsf orms • str earn 
functions of incompressible flows into stream functions. r.o.*f . 
compressible flows. The 'formula suggests the - possibility . of 
carrying over various physical laws which govern the motion 
of an incompressible, fluid to the case of. a compressible 
fluid. ■'••• 

.In a companion paper this formula will be used for. Con- 
structing a subsonic flow around a curve which approximates 
the boundary of an obstacle given in the xy-plane. (See 
NACA TN No . 973 . ) 

Another application Of the above result is to "distor- 
tion theory".! ' that is, the study of how the properties and 

1 The p os sibili t y for generali zation of the formula f or 
the case of • a mixed ( i . e . , ’part ially subsonic and partially 
supersonic) flow is 'discussed in- the paper. It is observed 
that for M< 1, £ is a complex- number , ; f or M > 1 a. 

purely imaginary number. Ther e‘f or e , f or -M< 1, f(£) ,1s 

a function of one complex ’.variable , tfhile for. M >.l it is 
a function of a real variable. 
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the shape of the. boundary change (in applying.' the preceding 
procedure, retain’in both formulas (4) and (6), the same 
function f) upon passing from a flow of an incompressible 
fluid to the corresponding subsonic flow of a compressible 
flp.id or upon changing the density-nre ssur e relation pf the 
fluid.' 


I. THE HODOGRAPH METHOD IN THE CASE OP AN INCOMPRESSIBLE PLUID 
■ i .■'■ A ■ General Representation for the Stream Function of 
Plows of an Incompressible Pluld in Terms of an 
Analytic Function 1 of a Complex Variable 


A stream function of a flow of an incompressible perfect 
.fluid .is a :harmonic function - that is, a function which sat- 
isfies' 'the' .Laplace equation . ... 

■■ + = o ■■ ( 7 y 

dx s 3y e 

Conversely, a function which satisfies equation (7) can 

be ..int erpr et ed as the stream funotion of a suitable flow. 
Since the imaginary part of an analytic function of a complex 
Variable satisfies (7), and for every function satisfying (7) 
there exists a function f(z) such that 


tyCx.y) 


j 


la f ( z) 


( 8 ) 


iB) is the "general formula” for the stream 'functions of a 
flow of. an incompressible fluid,. Here f(z) -. ranges over 
the totality of'' analytic functions. 


In connection with various problems in fluid dynamics as, 
for example, jet problems, another method of attack was 

" l In many instances an analytic function of a complex 
variable consists of several (or infinitely many) branches, 
each of which is defined in the whole xy-plane. These 
branches cover the plane many times,- Since a flow covers the 
plane or a part of it only once, each branch gives rise to a 
ph'ysically possible stream function. . However, here and' in 
th'fe "following* a'- function is always spoken* of rather thana 
particular branch of it. "' r ! ' 
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deveh-bped. the ba's-i-c" idea- of : which is t a consider' the flow not 
in the physical plane but’ r in the hodo graph -or so-called' loga- 
rithmic plane 1 - that -is t o-' i'ht'r b'dtice' as independent 7aria- 
ab r le"s" the component a' - Vj , v s bf-'-'tKe velocl-ty vect or- * 

y’ = ve^;- b v - + iv e = "5x . Wl =■ .<J>.- + 

and iog' v and. "" 0 , r e spect'ively , in§t ea'fL of . x and . y. 

This approach, leads to another .general formula whic.hr 
while it is more compli cat ed than - (13), has the advantage of 
being capable of generalization to the case of a compressible 
fluid. '• 


t i on 


In the case of an incompressible fluid the stream func- 


\j/(v,e) = \y x{iog v,e ) , y(iog v,e) 


• ($) 


is again a harmonic function of log v and 0 and there- 
fore 1 


\j/( v, 0 ) 


- 


Im f ( log v 



do) 


yields a ’’general formula" for the stream function (consid- 
ered as a function of log'v and 0). The representation, 

Mr( log v,Q) = constant, for the streamlines (in the logarith- 
mic plane) of the flow is obtained immediately from (10) . 


2. Passage from the Logarithmic Plane t o,, t tfee. Physical Plane 

The fact that the flow is considered In the logarithmic 

plane instead "of the original physical plane introduces 

" \ ■ • ■ ■ 

the transformation' Z* = log Z the author passes 
from the. holograph tg.th.e logarithmic plane. In. the follow- 
ing, in'many instances, it. is necessary to pass from the 
holograph to the logarithmic plane anl vice versa, often' 
without explicitly mentioning it. .This fact is stressel here 
in order to avoid coifusidn. ■ Thgt plane the .Cartesian coordi- 
nates of which are ..log v and 0 is denoted as the logarith- 
mic plane. 

See Hot at. ion, remark 1 . .. , 



14 


NACA TN No. 972 


dertain complications. In order to overcome them it is 
necessary to investigate more thoroughly t he - r el at i ons which 
exist between the flow around a given obstacle in the phys- 
ical plane and its image in-.the hodograph and logarithmic 
planes . 


Suppose' that the stream function ij/ = \J/(log v,8) of a 
flow (of an incompressible fluid) in the logarithmic plane is 
given. The following procedure yields the streamlines of the 
corre sponding flow in the physical plane. Since 

y = ilK ( see reference 6, p. 32), it follows inversely that 
dz 



( 11 ) 


Writing r. ve“^ k and noticing that the integration occurs 
along a streamline, \j/ = constant, and therefore d\|/ = 0 
gives (ll) written in the form 



( 12 ) 


Using the relations vcj)^. = -\j/g and (f>g = vvj/^ and noticing 
that along a streamline d\|/ = \j/ v dv + \|/g d6 = 0 and there- 
fore) d8/dv = ~\j/ v /\)/g. gives 


x + iy = s 





e- 


dv 


By separating the real and the "imaginary .part s there 
tained a. paramet ric representation, 



i s • ob- 


( 13 ) 
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for the streamlines in -the physical plane. One of the 
streamlines or a suitable part of it can be used.. as the 
boundary of the obstacle. ■ 

By employing the hodograph- method in airfoil theory it 
is convenient toiuse.;the fact that' the approximate form of 
the hodographs of the flows around airfoils of certain shapes 
is known. If It is assumed that the domain H which repre- 
sents the- ima-ge of the. flow in the logarithmic plane is given 
(see, e.g., fig. 2b), it is po s s-ib.le. t.o" construct- at first.. . 
the harmonic functi;oh: ,;-ii/(log -v , @ ) which assumes a constant ■ 
value, on the boundary h of the hodograph' ! and 'has .'the pre-t. • 
scribed' behavior' a ; t-:the point which corresponds to z = po. 

Then with (13) the form of the airfoil in 'th^ - physical plane 
can be' 'det er mi ne'd .■ ; . 

As is well known, for a given obstacle and a given angle 
of attack there exists a whole family of flows." If the ob- 
stacle 'has a sharp edge; as occurs in the case of an airfoil, 
all solutions but one have an inf init e . velocity' at the sharp 
edge. The Joukowski hypothesis consists ’ of the assumpt i on , 
that this exceptional solution, which has an everywhere finite 
speed, represents .that .flaw which has physical significance 1 . 

The- hodographs of the flbws around the same obstacle 
(in the -physical plane) lead, in general, - to quite different 
pictures: ’in , the hodograph ,-. and in the logarithmic planes. for 
ih'stance, in figures la, .lb,, and ‘2a ; 2h,- two 'different : flo,ws 
ar ound -similar airfoils .are. dndi cat's d . As mentioned, before, 
the hodograph. •' of .: Joukowski flows ■ lists', in general, a shape . '] 

similar to that indicated £n ■ f igure'^b (It is . noted that ' 

this domain is partially twice covered.) 

- ' " i ; ’ t _ . f 

If the hodograph method is used .i.e.; pbta'i n the Joukowski 
flow around some profile, it is at firs-t necessary to deter- 
mine the function 4f(log v, fi) which is defined in the domain 
E, and has a & .o u b d e t at A, the point which is the image 
of z = co. In order to- construct the stream function 
\|/(log v, 8) , proceed in the follo'wdng way.-: -Determine the 
stream function in the upper half plane '( Z-plane ) - that is, 
a function g(Z) which assumes constant values along the 
real axis and has a combined vortex and doublet at some point, 
say at '"Z ='i, and . then -the - fu&ct.i on . Z = Z(log V), which 
maps tbe Upper "plane into the domain ' 'E ’ traftsf orming Z ■ = i 
into the branch' point >• A of S . As.Wiii be sSen, a : family, 
of 'solutions is obtained for; this ppobifem. ■’ ’ - " 7 t 

1 A flow fulfilling the conditions: o : f. ihe ; v Joukowski, ,V! : ; 
hypothesis is termed a “Joukowski flow.” 
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The axis of the doublet Is assumed to form the angle a 
with the real axis to obtain for the complex potential with 
the circulation r the formula 


w(Z) = m log Z-zJL + _*L_ e ia - —Si— e" ia (14> 
2 tt Z + i Z - i ■ Z + i 

; : s . 

Reirfqirk t" ; The.term ‘i£ log yields a purely circulatory 

' ’ ■ 2tt Z + i 

flow' ('see 'fig'. 3, a’lso reference 7, p,-326); ifheredd 

mi L_ - represents a flow 'without any circulation 

\Z - i 2 + 1/ ; 

and with a doublet at Z = i the axis of which intersects 
the positive Z-axis with the angle a. (See fig. 4, also 
reference 7, p. 202.) 


The question of how to determine the mapping function 
has a more technical character and will be considered in the 
next section. 


• '' Suppose now a function \j/ =»‘Im£w fz(V)^] which assumes 

a constant value on the boundary _ h or H. (H is the image 
of D in' ‘‘the hodograph plane,). "TLh© boundary curve of the 
obstacle' "i s obtained if, starting "from sqm. point, say B of 


h , x' 
ever, 
curve 


and , y ’ are determined by, integrating, al-ong ■ h. How- 
in. general, the obtained curVe.' will, not- be. a closed 
, ; .In" order that this be so, '..it, is required that 


/ 

hi 


cos 


Js C 


+ V 


;a. 


t c 


4 


2“i ; 1 r. ■: 

2 . 


= 0 


(15) 


f 


h' 


sin .0 [\|/q' 


.» 2 
+ .. V . 


1 




dv’. 's 0 


(16) 


where'ithe integration '.i s carried, otit _ aloh^ the- boundary cuhve 
h < o-f = 4 H-J -Thus it is seen that , 1 n* : order that the obtained 
boundary in the physical plane- be- a ’ c’l osed purve, it ‘is neces- 


sary to choose 
tions (15), and 


r. m, 


fehd' : 


-a 


( 16 ) * arfe • sat isf led'; 


in such'"'a way' that both equa- 
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Remark : -In connection, with further applications for a com- 

pressible fluid two separate expressions have been derived, 


one for 


■the, other for y. Clearly, in the 'case under 


consideration they can be combined together, and (15) and (16) 
are then equivalent to ' 


*dwfv) 


’ dw(Z) 


( Z )dZ 


dZ 


(17) 


Since w(z) and V(Z) are analytic functions of a complex' 
variable Z, which are . regular in the- upper . half plane ex- 


cept at Z = i, and 


■tU) 


does not vanish there, (17) 


equals the residue of L^jL at the point Z = i .< 

. (Z) 6.Z ^ 

Write Z x = Z - i to. obtain for w and v, the series devel- 
opment S' 

w = 1 X- log S + — — + . . . , V = v a .+ T Z x + ... . . (18) 

Z, , 


and therefore 


■[- 


iZ| < p, p sufficiently small. 


" mie ia , ir 1 , - . 

— — ; ^ — T“ • • • 

Zi 2 tt Z ^ 

T" az, ’ + v x z 1+ y 3 z* + . . . 


1 dw 


, ia 

mie + IT 


Zl T « ' 2tt Z^ 
ia 


z, . JLv^ 1 j 


= „ + (_in + 


f. 

Y 


'CO • oo 

,ic 


l mie ia ^ i . 
” ^ * 


(19) 


Thus the-. above’ condition becomes 


m, tv i ct» 

T • Y x me ^ - . 

— + — — * at 0 

2u 


( 2.0) 
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3. The Determination of the Junction- Z = Z(V) 


Which Haps the linage of the Flow in the 
Hodograph Plane into the Upper Half Plane 

If the domain E is prescribed, then the function 

3 (7) which maps E can he obtained using one of the known 
methods in the theory of conformal mapping. For instance, 
•Theodor sen *s method (see reference 8) may be used to deter- 
mine the function which maps the circle into H and then 
compute the inverse function. The theory of orthogonal 
functions also yields (see reference 9, chs . 71 to IX) a 

i .■» * * ■■■ > 

simple formula for the function Z(V). 

f . Denote by c?y(7) a complete set of orthogonal func- 
tions. Such a set can be obta ined, f or instance, by orthog— 

onalizing the functions C <J(Y — a) ) where a is the 


■ u /- f l 

onalizing the functions 1 \ v(V — a)/ f, where a is the 

branch' point of the domain H. 

00 ___ _ 

3y K(V, T) = f Py(V) <p v (T) is denoted the "kernel 


function" of the domain. Then the function which maps the 
domain H into the unit circle, mapping the point A on 


the or i gin , is </ rr 


K : *(V, A) 

An (A, A) 


and therefore 


Z(7) =•- i 


•’/rTKjCv'. A ) - /k h (a , X) 

./WK*(T, a)+/k f (a,a) 


is the" required function." 

Remark: Equation .(10)’ may be written in a little different 

f orm . Writ ing '■ 

Z = log v + i (n — 9 ) 

gives 

^ (log v,9) = Im g (Z ) 


where g(Z) = f (Z — in). The passage from 8 to n — 6 
means that in the hodograph plane' the domain with respect to 
the imaginary axis is reflected. 
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’ll. THE holograph method ih the case op 

a COMP.RESS-.IBLE PIU1D - SUBSOHIC CASE 

■, ; 4. Introduction' 

V ■ * , • s * . * / 

' In tills part the h-o'dograph method' will he generalized 
to the case of a compr es s ible' f luid. 

The stream function \J> (x, y), in this case, satisfies 
a complicated nonlinear partial differential equation, 

(See (26).) If it is assumed, that the density pressure re- 
lation is. of the fnrm 1 

‘ P = P('P-) . 

where p(p) is a function of p alone', then the use of 

\|/(v,9) instead of • ^ ( x , y) . ( as-Chaplygih , reference 1, and 
Molienbroek have 'shown) represents an important simplifica— 

. tion. If the variables v and S are introduced instead 

*f x and y, the function ^ satisfies a 1 inear partial 

differential 'equation S(\t) .= 0. (See . (30) instead of ,a non- 
linear one, (26).)* 

Remark : It will he assumed that (unless the contrary is ex- 

plicitly stated) 

p( p> = A + c'p k • ( 22) 

where A, cr , and k are constants. However, the method 
developed here can he employed in the case of a much more 
general pressure density relation.- 

In the case, of an incompressible fluid, ins t ead of 
merely a statement that the stream function \J/( log v, 9) 
satisfies the Laplace equation, the" general formula (10) 
was given for solutions of the Laplace equation in terms of 
an arbitrary function f ... of one variable £ = log v— i9. 

'a . • * 

The main purpose of the second part of this paper will 
he to give an analogous formula f or a compressible fluid, 

■* and to derive from it the representation for' the stream 

function in the -physical plane. •. ■ 

j> As will be proved elsewhere, this result leads to a 

construction of a- flow around an obstacle approximating the 
given- obstacle (in the phys 1-cal - plane ) , - , .„ 
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After a short dris cus s.ion in section 6 about different 
types of differential equations and development in section 
7 of properties of the nuxil iary . f unct ion \(v) which is 
needed in the following, an operator is defined in sec- 
tions 8 to 12 (see (55,)). which transforms functions f (s ) 
of one complex variable s = \(v) + i0 into solutions 
\|/(v,0) of S 0 (>M = 0. Squation S o (\|0 = 0 is the equa- 
tion for the stream f unct ion . (in ..an appropriate plane)' 
for the case of compressible subsonic- motion. Then, if 
the following formulas 'are used 


X = X (v ) ~ 


y =■ y(v) = 



Po cos. 0 [(!.- M 2 )vf/ e 2 4 v 2 ^ 2 ] 
, P dv 

pv 2 *e ' 

pqsln 9 [ (1 - M g )4<9 3 + v^ y 3 ] 

• pv£ 


> ( 23 ) 




which are derived in section 14 and represent a general— » * 
ization of. (15). a parametric representation for the 
streamlines in the physical plane of the corresponding 
flow is obtained. 


In section 14 are determined the conditions that the 
image of the. given hodograph yield a flow in the physical 
plane around a closed curve. 


5. Differential Equations for the Potential 
and Stream Functions 

From the continuity and irrotSt i onal ity of the motion 
it follows that for every flow there exist two functions, 

<t> and \|/, the potential and stream functions, such that 


Po 


■dr 

■y 


X 


= V 


~P 0 




■\J/ s: <}> 

^x ^ y 


= V ■ 


(24) 


Here 'Vj" and v E are the Cartesian components of the 
velocity vector and-, p isHhe density. (See reference 6, 
pp . 228-^2,9 ''or- reference 2 p . 2.) From the. Bernoulli 

relation. l/2 v s + I ’= -constant , where . 1/ g = kep^ V g (k— 1 ) 
denotes the pressure head (see reference 6, formulas (13), 
p. 215, and (10), p. 214)., it f .oil owe, that 


1/2 V s + ierk (*. — 1 r 1 p k ~ l ■*:crk(k 1 ) 


\*r+ X 


Po 


k-x 


or ■ 
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V , P- 

: i 


= p ,1 r- 4 (k-l)k X ct i p I k v £ 
1 <3- r i a < 


i/k_: 




= P, 


= P, 


[* - i a;- s (^-i>v 

t 1 - i 




i /k-i 


y ( 25 ) 


a” S (k-l)(4> x S + 4> y s ) 


/k-i 


. v- ■ 

If (25) is substituted into (24) there is obtained for Q and 
\J/ & system of two nonlinear differential equations 


\J/ = $ 

y x 


O.-i 


a~ S (k-l) (t> 3 + <i> 3 > 


+ x = “ * y j>.l- I »o*0=-l>W** + < V> 


i/k-i 
i /k— i 


V (26) 


J 


It is noted that in the case where the motion repre- 
sents an adiabatic process, A = 0,. and in the case -of air 
k = 1.4. 

Eliminating- 1 ;. 4> gives for • . 



> (27) 


Similarly, eliminating . from (2$) gives f or- „ 4> : 


r~ - 

aS [$xx + ^yyj ='$ 


a 3 « a 0 3 


♦»X 'I-*?-*,, 

J r- ’* ■, .r ft f 1 ► . -f .■* ' 

r. ... J% *. . *. r • - L • c •. * ■ 


2 

X y XI 


~ ( k - l)($ x 2 + 0 y s ) 


I* 4 


(28) 


- • •' 1 33i?e .d'erdVat iom -of,, (27,):'. and .( 28<) i»- o$.l.t t ; erdv.her.:e . 1 ” The 
iequati'oa- (-28 ) ie, der ived in r ef ervemie; . 6 ., . p ■, \Bj30> s ;• .Conperning 
( 27 ),,.-. s;e© reference 10,, p. 5. ;jw'‘ '■ ; 
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As Chaplygin and Molenbroek showed, if the variables 
log v and 0 are introduced instead of x an-d y r Ahen the 
equations relating the stream and the potential functions 
become linear. If 1 



$ (v, 

e) 

= <f> [x( V , 

e), 

y(v, 

6)3 

^(v, 

e ) 

> 

1 ) 

V* 

II 

e)r 

y( v 

o)3 


( 29 ) 


is written, then, instead 2 of (26), it follows that 


f -Hiog t) “ °' ( t ) U " ^ + 


> ' 


&\|/ 




( log v) S( log v) dv 


> (30) 




Here 


•'K * v/[a 0 s - i (k - l)v £ ] 


i/o 


( 31) 


is the local Mach number. Eliminating _ <}> gives for \J/ a 
linear equation ‘ 


S 0 (^) 

. P> 



d\t 

5( log v ) 


}■ 


(32) 


In section 11 a general expression 3 vill be- given in the 
subsonic case (i.e., for M < l) for the solutions of (32) 
in terms* of an arbitrary function of one variable f. That 
is to say, an expression will be obtained, involving an arbi- 
trary function of one variable f such that for every f 
the obtained expression represents a solution of (32), and 
conver.sely every solution of (32) which is regular at the 
origin can be represented 'ih'.'the afore-rm.eu.ti on ed form with a 
suitably chosen f. 

^ee H o't at iony r e ra ark 1 . : \ 

2 A detailed derivation of (30) is given in sec. 3 
ref er ence 2,. 

3 1 '* * j 

It is noted. that, in-.-.s oc.. ..9 ; equat ion (32) is simpli- 

fied slightly by approximating the coefficients by ■ po lynom ials. 
The indicated result refers to this 's impl if ied equation. 
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Remark : It-, is noted that for k = ir Ab, K <1, oq.ua- 

t i oh. "(3'2)' hecomies- (i»< appr opri'at-e --variables ) the La-place 
equat i on . . 1 ■ 


. According to (25) and (31-)', 

P = p 0 Cl + (v/a 0 ) S 3‘ 1/S (33 ) 

M (t/* 0 ) [1 +. & ..or, );> ( T /a 0 )= M[l-H s r l/S (34) 

Since .1 - K® = (p Q /p) 2 = 1 + (v/a 0 ) 2 ' 

1 + (v/ a o) 


and 


where 


A* ' i- 


[1 + (v/a Q ) 2 ] 


\ . 

d ‘ 


d log v d\ 


X . 1 log L 1 t_ (*.l ! f ! - 1 
■ 2 


I7T 


Cl + (v/ a 0 ) 3 ] + 1 


equation (32) becomes 1 


W *xx * 0 


(35) 


..... 6 . A, Remark on D.iff.er.an-t.. Types of Equations 

The firsi purpose of th>e se bond" par t of this paper 
is t O;. give , a- formula for solniions of (3-2) in terms of an 
arbitrary function of one variable. 


Before the derivation of this formula is considered, 
it is well to 'discuss in some particularly simple cases 
the "general solutions" of this kind and indicate some 
characteristic features .of sujcb-.f ormulas . * 

The following three equations 2 will be considered, 
where, a and £ arp .real quantities: .... 

' ■ u. ■■■ — — i -| ! — ; ■! i _• 

x This result was f ir s’t'' ; "o'b ta ined by Chaplygin'. (See / 

reference 1, p. 99.) ‘ ’ 

2 * 

A differential equation Au^ + 2 Bu^.q + Cu gg + Du^, 

+ Bug +. P = 0 -is sa.id to be of;' elliptic' or hyperbolic type 
in the domain R, if AC — B ,fe > 0 or < 0 in R, respec- 
tively. .. . . * , . ; . . ■ - • ; - 
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JL (El + £3lf\ • = 
4 Vdp, 3 dd 8 / 


J>X = 0, 
at H 


= p + id , £ = 


M- 


& g \t/ _ 7 


dp. dd 


0 * •<>’-»> 8-=e-» 

In tho first case the "general solution" is given by 

^ = f ( t ) + g(£) s f(p + id) + s(p- — ia) 


id - •' (36) 
(37) 


(38) 


(39) 


in the second case by 


^ = f ( p, ) + g( d ) 


(40) 


where f and g are arbitrary- (sufficiently many times dif- 
ferentiable) functions of one variable. As p, and d are 
real variables , it is seen that in equation (36) there is an 
arbitrary function of one complex variable, and in equation 
(37) two arbitrary function's of a real variable. (Clearly, 
in equation (36) in order to obtain real solutions, for g 
must be ohosen the conjugate to f; thot is, f ( p, — id). 


A quite different situation 
•■■'tion (38). ...By t^e transformation 

reduced to the form ~ + — 

dd 


this cas e i-$ 


is 

\ 


0. 


me t in t he case of equa— 
= J 1 — p, , (38) can be 

The general solution in 


^ = f( i /l,' - p.' — d) ■ + ' g( i •J.x'Zr p, +■ d ) (41) 


It is seen that the arguments i */l — : p, ± d are complex 
for p, < 1 and are real for . p, . 1 . 

•The solutions behave quite ■ differ ently than in the pre- 
vious cases. It may happen that a. solution which is real 
-for p, < 1 becomes imaginary for > p, >1. Consider, for in- 
stance',' the- function 



On the other hand, there also exist solutions which renUin. 
real in tho whole plane — that is , 
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On the other hand, there also exist solutions which remain 
real in the whole plane — eiaople, 

( i (i - n ■+’ (id- p. ) +.*/)• .i ,2../(ji~Y ”i '+ f 2 ) 

• „ , Equat ion : ( .32 ) is, of mixed type 1 and therefore a situa- 
tion ekists of the type exhibited in ( 38 ) . Of course, the 
behavior is more complex than in the latter case, because 
(32) is not the simplest case of equations of this type. 
First, the function X(M) must be det ermine*!} which- may be 
done by reducing the equation (32) to the canonical form. 

7. The Function X(M) ' 


In this . se.c.t i.bn the. function X(M) is intro- 
duced. For convenience , an intermediary variable H = K(v) 
given by 

dH ( 11 = JE. (42) 

.... dv v 

is employed . 


The equation (32). becomes 

= ("-Y (l-M*) I |-^f] ■ + |~ = J (H) + Y! = o (43) 

\p/ Ls9 3 J. : Sh 2 de 2 8 h 2 


sW 


If M < 1 - that is, in the case of subsonic motion - then the 
coefficients of both ty’gjj and are positive and there- 

fore the equation is of elliptic type. If M > 1- that is, 
in the case of supersonic mbtion- the f or ego ing , coef f icient s 
have different signs, and', the equation is hyperbolic.: 

i * . • 

In order to obtain’ X(ld) , (43) is reduced to the so- 

called canonical form. (See reference 11, ch . I, sec. 1.) 


Intr oducing 


i = A (H) + i0, i = A (H) - 19, 2 (44) 

where ’ ' 


dA ’ 

dH P ( 


d -1 Vl - M S = */ 1( E) •, that is 


dJkl y) J 


dv 


=' v 1 (l-li 3 j ' l/2 ( 45 ) 


!The local Mach' number M plays,' In tFe case of ("32) j a' 
role similar to that’ of jj, In the . case . of. ( 3.8 )-.. ,. • 

.. ®It' is not ed 'that., for.. M < 1 , ,£ and, T ..are complex 

quantifies ' whieh.are conjugate to each other, for M > 1 they 
become two (independent) purely imaginary quantities. (See 
(48), (49).) 
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the equation (43) becomes 

l„w , , 4 !^ + i_£i_ v||_ jj 


(1 - M 2 ) 


l/2-j 


-..,4. 1 - M S p 


JLat HJ 


where 1 

N = N( £ + I) = N(v) 


= + m + n f|£ * ?Al 


(46) 


(k + 1) 


v 


‘h s - 1 


(k + 1) v J 


-3 /a 


- -(k - 1) 


(k + 1) 


M 


8 (1 - M 3 ) 


^-x/a 

J 

(47) 


The function A(v) = ^(M) may now be easily evaluated. From 

■ji/a 


d A(v ) 
dv 


' a ’ I 

= 1 L 0 w 8 


( k + 1 ) v ; 


-(k ~ l) y 


f/a- 


i ' 


v iv - ~ 

L 0 2 

it follows, by a purely formal c-cmput at i on , that 


A (v) = A(M) = ~ log 


2 1 / 2 

1- (l-M 3 ) • 

/ E 1 Is X 1 

/■l +■ h(l ~ M ) \ 

_;i + ,(i m 2 ) 1 / 3 

\ 1 - h ( 1 - M 2 ) 1 / 3 / 


(48) 


h = 

. \ k : + 1 / 


- iV^ 3 

' for k > 1 


For M < 1, \(M) is a real quantity, if M > 1 , t hen 

' I * *. , ' 

A(M) = - i 


(tan -1 ) Vm 2 1 -■ i(tan -1 ) (h ~ 1 )J / (49) 


is a purely imaginary Quantity,. 
. 1 See Hot at i on . r emar k 1 !! 


For k - 1, . -X(.M) = — 1 o g 

3 '* 


1 +• (1 M*) ' J ■ 
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> •* » - 


£.S&aEJsJ .If k < 1, ■ and M < 1, 


A 00 = 7 (tazTV- L'yT^(i.-:M 3 ^ /2 ] + 


rsa/s 


i 0 g / 

f l Cl + ■ (1-H») l/ ‘V 


i“ * » 




y (50 


■7 




If k = — 1 , M < 1 


1 (H) = 1/2 lo:g f — 4~^~STT75 1 *'i/* 

' U+ J . l[l+ (v/a Q ) 2 ] l/s + 1 J 


or 


A 

M = ae^/Cl+e 8 ^) and (v/a 0 ) = 2e K /(l-eS x ) 


For the application of this theory, the inverse func- 
tion M = M (A ) often is needed. This can he determined 
either hy preparing once and for all a diagram A = A(M) 
for a fixed value of 1 k or analytically, representing 
M = M(A) in the form of an infinite series'. For M < 1, 
k>l there is obtained: 

• ; '> ‘ 

' 1 / * 

1(2 A.) = /l -” M 2 ( 2 A ) = 1 - X - l/2 ( 2k + l)-X® 

(4k 2 + 6k’+ 3) X-- l/24 (24k 3 + 6Qk 2 + 76k. + 29 >X 4 . 

-~(4Bk 4 ■+ 212k 3 + 392k 2 + S38k + 103)X 5 : . 

— — i-(480k 6 + 2976k 4 ' + 7968k 3 + 10788k 2 ' +' 7266k- + .1935)X® 
480 

-|^^(2 88 0k 6 + 23472k 5 + 84232k 4 +■, .162124k 3 + 173940k 2 

+ 98086k + 22675 )X 7 — ... ( ( 52 ) 


X = 2 


■ /TkTT Ti ' 

( <.*+ - u- i) i/£ Y (k - 1:) „. 

\(k + 1) 1/2 + (k - 1 ) l/s V 


If M > varies between 0 and 1, A varies between 
—oo and Oi The proof of the convergence -of .( 52) for A<0 
is ' given in section 15. ... ““ 


l The .corresponding values of- 2A, M, and' v/a 0 for 
k = —0.5 and k = 1-.-4 are given in the . tables • la and lb. 
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8. Tho G-eneral H epr e s enta t i on for the Solutions of (32) 


lhfi.-2r.aa (53). L et the f unct i on E(H, 0 , t), — 1 = t ^ 1 
he a solution of equati on 


0(e) = 



l H 


X 


A H 


+ 13 


+ a mEV + -MsL. . o 1 

2 Ajj "" \/l — t s 


(53) 


which has t h e pr££^r tjr that 


("/ S,.t 

lGi + 


A h Vi - t* 3 Jl - b Z A 


HH 


(54) 


coa t inuous at t = 0, a jfc 
S is. g ive n by (43 ) , s.n£ 


t (A + i0 ) 2t (A + 10) 

A = 0, and a.t 0 = 0. H ere 


A h S' A„(H) = */l(Fj 


l H 


(See (43) and (45).) Then 


'Kk 


/*i 4“ ^ r “| 

I,9)=y E(E, 0, t)f |^i.(A(H) + ie) (1 — t 3 ) dt/-/l - t 3 (55) 

where f ( s ) is a n arb itr ary , twice— dif fere nt iahle f unc t i on 

■21 .ana yahialiAa , it All In a s nlnAAsn °£ S(ij/) = 0. 

Pr 0 of : It is noticed that 

f H = p A H ( 1 ~ t 2 ) f ’ , f fl = ~ (1 - t S )f < , f t =-f <(A+ i0)t (56) 


2 

wher o 

Ther of or e 

- . _ I A H (l-t e ) 

E 2 t (A + 18) * 




f - _ i Cl - t a ) _ ifH 

e ~ 2 t(A+ ie) tf e " A H 


(57) 


How, "by 3 (55) 

x [ ]^. indicates differentiation with respect to t. 

3 

In order to he completely rigorous, the integration 
is carried out along the curve — l<t<— e, t= Ce lc P, 

— tt ^ cp ^ 2 tt , € < t < 1. Then the integrands in (58) and 

(59) remain continuous along the path of integration. By 
use of the property that (54) is continuous, it is possible 
to let e subsequently approach 0, 
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H 


-r: 


E H f 


dt 


♦r 


■£ v ^ 


n/l - t 2 

dt y >+1 E_ 


Ef g 


dt 


Jl- - t E 


> (58) 


1 7 1 ts A B 
t(A + 16) dt 


dt 


Jl - t 2 

Integrating "by parts giv es- -for - the last term in (58): 


J 


>+ 1 


+ i 

( E a/i- t 2 Ah VI 

kt -1 

f*E v'l.— t 2 Ag ^ 1 \ 

^•1-t 2 

t( A +’*i8)V 


LS'jtrf-A+ie) t __ 1 / 


(59) 


% 


In an analogous manner there is obtained 

+ 1 f l __±L_ , + f£ 

I /l i P l rt i. / A i .* A \ / 

t 


= r iVI 

J-i L Ji - t E \2 t(A+ 10) / t J 


it 1 g ■ dt—J.1 1 *■ 


2 t( A + 10) 


t = — X 


(60) 


S 


Now, differentiate ( 59)- with respect to H, (60) /with respect 
to 6, and multiply by Ag 2 . There is obtained finally 


\l/ + A \i/ 

V HH h ^ee 


-r: i 



*(! 

_ j” Ji - t 2 


' e H^H + ' 
Vi-t 2 

yipr) 



■ t 3 

L t ( A 

+ 10) 

E Vi - 

t 2 \ 

2 ( t 

K 

A H 2 fei 

/ A . A \ f 0 




s / eh . a hh e . i a iY| \ 

vi^ aft + — ; j dt 


A E f H 


+ K 8 f e A H 


2t(A + 18). 


EAh : 


[ s H A H - 

MZ. 


EAg* 


( A + is ) 


+ ^ A gg + EgAg - i 


(A /+ 18 ) 

| ^ Z A S ? , [a f + iA *f Q 
d.2t(A •+ iS) L H H H 6 






( 61 ) 
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Using tho.last relation of (57) it is seen that the second 
and the', f our th term in the second integral of (Si) cancel 
each other, and if fg = Ifg/Aji is substituted, the last 
term of (61.),. vanishes. Employing (57) again and integrat- 
ing- by parts' gives 



Eg yi'-T^Ag 
2t(A + .18) 


f t dt 




/Eg - t E A H 


\ 2t(A + 19) 


+ i 

f + 

— i 


f +1 . t at 

jL x \ 2t(A +:.ie) . /, 


(62) 


and 
+ l 


x; 


E 8 f 8% ' 


dt = 


-r 


+ 1 ' Eg iA H 3 t 8 ' 
. 2t( A + 19) 




lEeAaVi t» f | +1 + t dt (63) 

2t( A + i8 ) _! J_x \ 2t(A+i0) y t 


Using (62) and (63) gives 


\1/ + A 8 ^ 

¥ HH H 0 6 


= / + V 

1 V. 


SU.I 


vr~ t® 


A a^ EH + 4v + A HH E 

i/T=T» ’ L t( A + i9 ) H W 6 A n 2 


)]J 


dt + 


f A - t 2 
\t( A + 19) L 


Efl-Ag + Eg iA.jj 2 + EA 
» 2 


HH 


} 


t = i 

t = — i 


t- 

( 64 ) 


which implies theorem (53); 


9. A Simplification of the Problem 

Following the present line of attack the next step is to 
investigate the solutions of equation (53) arid -to. determine" 

those among them which are most appropriate for 'the develop- 
ment of the theory. 

However, the mathemat i.ca.l, amalys is of thid qllfest ionvhas . 
not yet been developed to t*he e-xlt.ent needed iff -the cas.e under 
consideration — that is, im- the'- : gene'ral case of &n equation 
of mixed type — and to work out this mathematical theory* 

1 The author will develop this approach in a future paper. 
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here would 1 ead' 'ou't s id e the s'cdp'e of the* present paper. 
Instead of this, two simplifications are made, "by which it 
is possible .to employ already known mathematical results. 

j’irst , only the subsonic case will be considered. This 
means that the solutions of the equation (30) will be con- 
sidered. only in the domains where the equation is of elliptic 
type. Secondly, function J in .( 70) will be replace.d by a 
polynomial in .which vanishes at X.=— co . 

•' In the. ‘‘pas e of an incompressible fluid where ^ is a 
s o'lii’t i.on 'of the ’Laplace equation there is obtained for the 
stream function thd representation 

^ ( v, 8) = Im f ( s ) , s = log v — i6 (65) 


in -terms of an arbitrary function f , of one variable. 
(See equation (10).) 


Generalizing this result, it is found in. .the following 
that the stream function of a subsonic flow of a compressible 
fluid, which is a solution of (32), can be represented in 
the form s - 


'K v , 8) = li® Im 
m — 

CO 


^h(v) £f(': 


a) 

ft- 


l ^ f : ,]}< 


n= i . 

where £ = \(v) + i&, -and ' -H(v ) and 
are functions whi-.ch d'ep end ' upon m, v 




■,Q ( m n) (v).,' 'a ‘ 2, 

and k. For 


k = 1.4, q = lim- Qrf . are -graphically r epr es'ent eu. ; in 

.% m- '=■'*' • - " 

table-- lb . 


The remainder of this section and sections 10, 11, and 
14 ar e • d.ey=o-t ed , .t o , an exact- -formu-lat ion and der i-vp. t ion- of the 
f or ego.ing^ rppr es.entat inn for the- stream- f unct ion’-. 1 11 ' ■ •* "■ 


'J -> ‘ v . -i : ■' 

Ip,, Pf dear t-o link the .ensuing 
elnatical procedure', equat ion, >( 3 Q ) 
form L 0 (ty) = 0, (see ( 46 ) > , by 


analysis - with standard’ ma'th— 
is reduced' -to - the canonical 
introducing the variables 


1 The mathematical details of the proofs in secs, 9 to 

11 may, for the most part, "be omitted by tHe reader whose 
primary interest lies in the field of physical applications. 
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i = A ( M ) + 10, £ = A ( M ) — 10 (67) 

» . r , « ' | ( ( , r * ' •- 1 f ■ / • • I 

(See (44) and ' ( 48 ) • )' The equkt ion , ( 53 ) * "becomes . 

> .* • t 

® X (B) = (1 - t 2 )(S^ t + NSj) - t~~ 1 ( E £ + HE) + 2it L 0 (E) =0 (68) 

( 5.4j) will te satisf ied 1 if ( E 7 + NE)/t£ is 
. po-.lnt t = 0 , t = 0 . ■ • . . . ’ • ' 

If now, instead of ^ and E, 

ty*'- = 4t(I+D-*|\fr and E* = Nd(C + £)}jE (69) 


She condit ion . 
regular a.t the 




are cons ider ed, then L 0 becomes 

. ' ' * L(\|/*) = + E^* = 0 

where ■ _ ; ■ 

E = - ( N £ + N 2 ) 


(70) 


( k + 

l)v 4 ] 

! _ y 

[ a ° 3 “ -jg- ( k + 1 ) v s 

1 3/2 | 

a 0 ® — ■— ( k — l)v 2 .. 
L & - 

l/2j 

f 2( 1-M 2 ) 1/8 

V 


(k + ' l')'y 8 


’ '64[a?~ ^(k + l)v 2 ] 3 [a 0 s - i-(k'-l)v 2 ] 

•- .( k . + 1 ) T^lfc a 0 4 ,'+ 4 ( 1 - 2k )a 0 a v g 4r ( k +_1 ) v*3 ' - 
64[a'o® - (k + ljv^ao 2 - -|- t (k-l)v'^‘ : '' 

( k + 1 ) M 4 [ - ( 3 k - 1) M 4 -4(3'- 3 k ) M 2 + 16] ’ *’ 

64( 1..- M 2 ) 3 


(71) 


Since A, and M, ar e: conn e:ct ed- by the; :relat ion (48), 
the expression E is a function of A which has a pole of 
the second order .at A' 0.- . And-. . A = 0 lie's on the boundary 
of the int erval- of:-, v-ar-iait ion of. ' A, since, ,if- K ranges over 
CO, 1) , A.- -ranges - -ov-.er •(•t-rep', Q)-. . 
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The second simplification is made by replacing the 
function by approximating! function . F which is a 

polynomial 'in e 3 * of the 'order m . and' vanishes at' X = -w. 

In -section 16" it is -pfoved 'that in 'every interval 
(- 00 , X'o ) , '-Xq < 0,-, trhe original F may he approximated 
arbitrarily closely by such a polynomial!. This means that 
to every X 0 ,< 0 and every e > 0, ? ther.e is de*aayftiflied. *a 
polynomial y m in e \ F m (~=>) = 0 -such that 


F(X) - F m ( e X) 


g e 


for _t00 % X X< 


(72) 


The following is now proved. 

~ L am‘i£a- - (67) . T o every polynomial :. F m (2X) in -there 

exi st s a constant c , such that 

’ '|d% m (2X)/dX K jSc(E+l)i/(-X^ f for X <0 and K'= 0, 1, 2, '.(73) 

Proof?' ' Si nine 


m 


’Pi 


( 2X ) = ^ c^e 2 ^^, Cg constant 


K=i 


it suffices to prove that the inequality (73) is valid for 
thfe derivatives of a single t erm-- ! -e s ' >v . But d K e s ^/dX K 

= s^e and as K — > a> % ! / (_he ) coj. therefore there 

exi sft-s v a - d'hbh'^that i ... - 


d K e S X / d'X K ‘ c(K+ !)'./(- X) K t 2 , K=0, . . . 


(74) 


1 1 n some instances it is expedient to approximate., F v 
by the, sum of a polynomial and a function which becomes in- 
finite as at X = 0 . 

X 2 - • ■ 

In this second case it is necessary to use the result 
of reference 3, sec. 4, instead of theorem (83). 

, ■ ?. y 

Note that in (4,2) of reference 3, Of 1 (k + l)(a - X) k+S 


should read - ^ k + , ^ > in (4,3) [T (n" + l)] 

(a - X ) k+3 

[T(n + l)^ - 1 ; and on line 11 of p. 279 lira c 

' a ; ir • 


should be 
(n) 


n 


03 


should be 
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In sections 10 . and 3,1 an integral .representation will 
"be derived in 't eras of apalyt Xc ' function s for.the solutions 
of Ii'jji , the equation ’ ; r‘ 4 Suiting from replacing E "by E m 

in £(vj/*) * 0 ; that is, L jn (\j/*) = \J \i*~ + = 0 . 

10 ,. lemma : 

Lemma ( 7 5 ,) : . IX "E * ( £ , L t ) Xs. a>, .solution -of eauatl on 


-1 


8f 


then 

i' 


■L 

w.h e-r e 

* ( 

variable 


' e*(£, £, t) f t (1 - t s ))at/(i - t s ) 


+. v*] = 

0 ( 75 ) 

= O' and 

0~» 

II 

0 

- tV /s 

( 76 ) 


Proof ; Differentiating with respect to £ gives 

1/ 2 ■ ! \ ' 


( 77 ) 


4-t ' 


j' ‘S* fT-dt’ / ( 1 - t ~) 


( 78 ) 


• -1 


Differentiating again with respect to £ gives 

.+ 1 


\I/* ~Z = 

¥ £ t , 


/* E*_ f dt/(l - t 3 ) l/S 

•-'-l 


+ / -E* 3 

-L r. 


f dt/(i - t 8 ) 1/8 ( 79 ) 


If it is noticed that 


. fg.« -|(l t S )f' and ••ft". -• - £tf.' 


where 


d f ( s ■) 
■ ds 


3 s 4d ^ 


there is obtained 


and therefore 


'f I (l : > ;; t 3 ) f 

£' 2. ft- 1 


( 80 ) 
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^ rbfj - 4 if 2 ^ L ■*! * t4t 

=1 | jT + 1 t pi - t a ) 1/a _* 


(1 - t 2 ) l/s 
- : s >* ,f 

2 £ t £ 


t=: 
t = 


. 3 £.f B rl dt,31 > 


(■by integration by parts). Substituting the obtained value 
into (77) gives 


E* f 


" :+: .£ +lf {[- 


2 \ 1 /2 


t 3 ) 


- E * 

att T . 


t=i 
t=— i 

♦ itt 

J t ( 


l + F m E* \ 
1 _ t ? ) l/2 -’ 


at •( 82) 


which implies lemma (75) t 


. 11. The Representation of the Stream Function in the 

Logarithmic Plahe for the Subsonic Qase 

Theorem ( 83 ) , Let F a ( 3.X ) ^be gn analvt ic . ftmct ion of 
a r eal var iable , X , def ined f or X < 0 , which posses s es tho 
property . that , . * . L . ,• - . 


d E F 


m 


d X 


< c(K+l)l 
(-X) E+B 


for X ~ 0 and t = 0,1,2,... (83 


where c .is. a suit ably chosen constant . 

Further , let Q^ n ^(2X), n = 1,2, . . ., denot e a s et 

of funct ions which are def inod . by the recurrence f ormula : 


(2n + l)^ n + '^ + .fti x \^,) + 4F m q/ n) = 0, l)= - 4r m ' 

n ) (a) = 0 , a< 0 

Finally , let g( £ ) bg &n analyt i c function r egular in a 
domain B which contains the origin, Then 


( 84 ) 
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r 


^ * ( X , S') = Im g( £ ) 


y q (n) ( 2 X) r . . . din. .dtj (3F> 

2® n n! o ••'o J 


n*=i 


will be a solution of 


± A\|/ * + !-( 2X )\|/* — 
4 


= 1 Q 2 ** ^ d 3 **') 


b\‘ 


be- 


J + F m (2 X)\|/* 


s\|/ * -• + F ffl (2X)\|/* ■=. 0 


( 86 ) 


which _i.§. def ined in every s imply connect ed domq. in lying in 
the int ers ect ion of H and 3 , wherj • H . denot es the domain 
0 3 < 3 X s , X < 0. 

.Ex.fi.g.f : If 

. .03 • ' 

. E* = 1 + t.t 1 Y (t t l/S ) ?ri '^ 1 Q (n) (2X) (87) 

, ■ ; ; : . n-'i. 

then it must Ihe shown that-.E* satisfies the equation 

G g (E*) = (1 - t S ) E± ^ - t _1 E* + 2 £ t £ E*^. + E m E*J = 0 (88) 

(see equation 76), and that E? / £ t is regular at £ = 0, 
t = 0. E* formally satisfies L tqhation (88). In fact, 


+ • . . + (n.-l) + , . 

b t ' 

s h * + i £ 4l } ) + • 


(89) 


Thus 
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-E* t = 




■- t 5 V 


< I » 


t s E* = . . . + (n~2)t sn - a t n ~ s q i [ n ~ 3) 

* ‘ , / ’ ’ *- . v ,» 

+ (n-l)t® n -i^-i Q ^n- 1 ) + . . . 

r * ■> 1 

•» . , ' V 

-i 

+ !*?■»-* + -• . 


r 


—2 1 ^ ^ 


. _( n -a) t aa - = . 




- - 


2 

-2t = -2tS3P m - . . , - 2t?*-Bta-iy ffl q(n-l) 

- - . . . 


or 

<> s (s*) = 


t rco 


■^x 


+ 4B 


m 


2a- 




A 


> (90) 






(n) • '. 


* *Sr 




(si) 


which, implies (84). Now, proceed to .the proof of the con- 
vergence of (87). If. A is a dominant of B' — that is, if 
for all derivatives ^ d^ifc/ d«X? i ' X = 0, 1, 2, . . . , it fol- 

lows that 

h ,-O.y A , 2 , . . (.-• “ a %. X <0 (92) 

' / , ■ * , ' ' " ' * *’ * ''•* ' 

This will "be indicated hy 

B <■< A. > ’ t»r. A >> 3 (93) 

By ’the r e:cur r-e el .ce„r :f o rmu la . . 


d 1 ^/ d£ E l^', mi d\ K i •• i ■* yX- 
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( 3n + l)Q ( x n+l) ?;Q ( V x ) .+ 4c.(^) 2 + (94) 

i , L*Ja‘ ’ ' J 

q(i)= 4c/(-X) s ., . , Q- n ^(a.) : ' - O' 


the functions 


(n) 


are iutroddced , 


Writing = c (n) (-X.r :(n+ ^ s ' 


gives 


(2n+l)c (n t l) (-\r (lx+ f ) ,>= ((n+l)c (a) +. (-\)” (n+2) (95) 


and from this is obtained 


lim 


n 


c ( n+i ) 1 

' J n) - •••"■ 2 


(96) 


Thus, the series 


1 + 'X t n q^ n ^(-2\) 


converges for 


— 2X 


n = i 


< 1 or X s + 6 s < 4X S 


(97) 


(98) 


It is shown now that 


q( 


(99) 


Clearly, the and all derivatives d^§'' 11 ^ /dX^ are 

positive, ..Further, -by ( 84) , • ( 94)’, and (83) it follows that 


. .$( i) <<q(-i)- ’(100) 

Equation (99) follows by induction. Suppose it holds for 
some n, say n = n ; then, by (84).'.and (94), 
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, u..jL1. _ _ x , 

. > v ‘fo^.+ 4 / dx 1 

( 2(j,+ l ) l x 4 -1 

^ — i-r < 101 > 

(2n+_l) L X J (-X) s J 

. Further , since, ; all derivatives of q”(|J' + l) are com— 

' ' | j 

binat ions of the derivatives q_ ^ of Qp ^ ' with positive 

' ‘ • ** i[i+i) 

coefficients, an.d| V pince the expressions for Q, ^ in terms 

(u,) „ - ' ‘ '■ .... . X.- 

of Q,^jr hay^the same str'uc/ture as th'o expressions for 

^XP +lJ in .totjaB x of , ( 1C °) follows . Eere Q^k =• 

'■* ■ •*- , y. ■ ’ ■ vV* *$' dX "■ 

*in (.87) is a solution 

af ( 88'} ‘s'a't'i-sfy ing* the' coh’dii ion : Ef/^t regular at t, = 0 


of 

and t = 0 




3y 1-emffl a ('7J?).- i,t>. Follows' that 
+1 

..... r. / 

'I'* = lm’ 




where ' f (.Z) : an arbitrary analytic Function'' of a- complex 

variable Z, is a solution of (.86) . The series (102.) con- 
verges : un if ormly for 0 21 '<5X 2 . Therefore af ter. replacing .. 
$*■ in (102) by the right-hand side .of (87);' the’ brier.' bf 
summation and integration in the’ resulting expression may' be 
changed to obtain 


*\ : 1 r> 


: * : r . ? { y .1 


\ t (i-t- & rj-<t7^;i-t- s ) 1/s 

00 ^ ‘ ■ * ■ * '* 

' [' + 1 n h2\)i n A t; 2tt f [i ( l-t 8 ) : dt /( i-t ^ i’/s\ ( 103 ) 

1 '' •"ifeY* ? ■ ' *JJLi "• .1*. ' ' : J • “ -*" ■ ■ v* ■. ■ v 


Let 




ay Z ' , lak'd ' wr it e, 1 


••■1 1 of 


n=o 


00 

1 r(-p) is the gamma function / 0 e -t X t-^ - 1 dt ; so, for 
integral valuos of p, T (p + l) = pi = p (p — l) . . 


. 1 . 
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e 


f f [-J-su-t 8 )] <it/d-t E ) = g[°ht) 

:«>(£> . V 

. J * 1 . . — 2 u (n+v) . . (v+l) r (u 


■i?=o 


n =.. 0 , 1,2, 


( v + l) 

(104) 


Then 


4- 1 

e [, 

-i 


t en f r i £ d-t 8 ) 

. U ^ 


n+u + 1 


y CUTV - y_l 

•Av- 1 / (l-t) t 

- s' , SU* J ° 


y j. X^M^s) 

7 .2^ . r(n+v+l) . 


dt 


v=o 


„ n+v 

t a» 


■ (-iX-t)- • ■ i'l r&k 


r ^ v+ 4 ) r Ci 


t 


n+v 


n+v ). . v . (\> + l)r(u+i) 2 V 

= l— 1 ^ ( 2n ~ g ) • - . - , 1 . g (n) (£) = _.£1.3n + l) -(n)(M ( 105 ) 

2 n S ^ 2 sa F(n + l) g . U; U05; 


Substitute the last term of (105) into ( 103 ) to attain the 
express ion- (85 )»•■•. „ ' • ' 

12'., The Evaluation of the Coefficients 

It was proved in section 9 that if F 1 is replaced by 
P m the series obtained for E* converges, and (85) rep- 
resents a solution of (86). 

\ , ■ ^ . < 

It is important for practical purposes to compute the 
Q,^ n ^ explicitly. 2 Since m can be chosen so large that in 


a In thi,s -a.nd t he .-f. oJllowing section advantage is taken 
of the remark in footnote ’2 on : p. 10.- 

3 

In a later paper the corrections to'. be made in order 
to pass from the obtained in the ab ove— de scr ibed 

manner to the Q, m ^ n/ obtained using a polynomial F will 
be determined. , It . will .be,,. se.en ..that., . in. general, these 
corrections may be neglected. 
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the given interval ^o) . X 0 < 0 , F m .and any required 

number of the derivatives of F m differs hy less than any 
prescribed- e >0 from the corresponding derivatives of F, 
th.e, will be - comput ed- us ing the. function ■■ F instead 


of F 
)( n) 


m ' 


As will be shown, the expressions obtained for 


Q' I1; consist of a finite number of rational and logarithmic 
t erms in . 


T = ( 1-M 2 ) 1/2 


(106) 


(For the relation between X and M see (48).)’ From the 
second equation of (84) 

X ’ ’ i(2\) 

d( 2X ) 


F - - - -- - - dT 

dT . 


v . ^ = •— 4 J F dX = —2 r 

'• * -i- £b " \ -<» . 1 .. 

= - if jiiitid _ ±2 e + s (?£~ 7 ) . + 4(k+3) 

L T S T 4 T s 


■ - ( 3k— 1 ) T 2 


J (T 2 -!) (l-h 2 OJ 2 ) 


dT 


= ( k+1) 
8 . 


- T 


3T~'" — k*l)T 


log 


„ fc-1 „ 

1 4 T 

k+1 4 


( k 2 — 1 )k 2 — 1 ~ “ 1- T 


'f 


. . : . . :fc+l : 

i — * ■ » ,• ■ , 


k— 1 
T 

T = 1 


(is?) 


Setting n = 1 in the first expression of (84) gives 


= -Q^ - 4FQ^ 1 ^ ' 

A. A. 


-f 


or 




X 

X 


•Xv 


J-sd 1 } (2X) , 


•' ; - 4 r fq,^ 1 > dX = ' 1 ) + f q (l) dq ( ^ 

^ — —03 OO ^ 


r - i + 4- of 1 ^ « 4 r + ~ Q ( 1 ^ 


q(?-) ; (_co) 


3 \ 


(io'8 : )- 


t v 




-J-’ 


1 . 
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It is noted that, in, general , i.t f qJ Lows- .f-pcm. ,£ 64 ) that 1 -- 


<2n+l)^ n+1 ^ =‘ 4 /%(«>, ,dA.;^ + r ■ dQ^ 

A J A J 

- Tf-») 

since from (d £ I'/dX K ) =0 it follows that 

\= — qo , 

; <>..> . 

^ n) (-«) * = °*' . 11 = °! 1* 2 > ‘ * * 


( 109 ) 


Thus there is otytE^ned 


5 q, 


1 n(0‘ 


( 3 ) = _ ip + - ~ f q' 

3 x .3 ^ • 3 J 18 


» —00 


or 


> - - + — r4 (l> - if f + (HO 


wber e 


re 

128 2 9 t 


48kT 2 


-4(3k-7)T -2(3k- 


:— 1 ) T^| 


and 

/ 

-CP 


„a _ 25 ( k+1 ) 10k( k+1 ) , 39k 2 —8 Ok— 1 15 

jj d.A *- + g 

9T 9 T 7 5T 

4( Ilk +48k 2 +59k+20). 41k 4 -fiSk 3 -^ k 2 +112k+91 


S( k+lOT 3 "' 


( k+ 1 ) T 


+ ii gk 4 +17k 3 -15k S -21k+8) T _ (3k-l) 2 ( k+1) T s 


(k-ii)-- ■ st- 


3 ( k— 1 ) 


In order to obtain from ( s e e ( 69 ),) It is necessary to 
have ’ , r ■* ■ 

.M 1 


H = exp 


( “ f I(4^)V- § 

^ j ' (l-M s ) lA L 2+( k— 1 )M S _ 


2( k— 1 ) 


( 111) 
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Al ( y , 0 .) = H [im eU) + y ^ 1 } ‘-y ^ K ^-.g (n ^n1 (112) 

• ! * ' l-.'L*.', 1 . 2 an r(n +i) J 


n= a 




where Im = imaginary part and. g^ n+1 '®(£) = f g ^(£)d£, g^ D 

, . r*‘ • * . . ' ^ 

- ' !l - • . ... . .. . 0 

= g( i) . In order to evaluate the terms £>f (136) and 

are needed. ‘ . v . 

Dif f er ent iat-ing (112) with respect to 6 and to v, 
respectively, gives ‘ ' ■ ‘ ■ : ■■ * f 


’ ; : * i ■ • r .. 

• = H'rite 


g 


■t * 

j (Re*,? ye&l 'part 


'+ ! Q. ( ^Re ; g + j 


A a) ( 

Q Re g 


l} + . . .] 


(113) 


\J/ T ■=■; R.( ° ) . . ,1m. g g + r( 1 ) Im g + jg-R^ 2 ^ itm g^ 1 ^ R 1 ' ^ Im g^ 2 ^ 


•r c . v . 


wher e 


; '• 1 ’ . + . y -\ R^ Im g^ n ^ + • • .(114) 

■ 2 2n “ 2 r(n) ■ . 

•: fi( °j =.*.H %r R- 1 ^ VH*;+ : .• . ■ 


dv 


d v . 




/„„(n ,-OV ,, 2n(2n-l) TT «( n )'l d\ _ _ „ , 

h — J . V ’ * ’ ’ 


The values of H-, Q v ■= • 1 im for k = -0.5 and'' those 


,(n) _ 

m— > oa 

of H, C„'‘ ru ,R vrtl1 * lim R m ^ n ?^'T,o.r *' 'k = 1.4 ■' are given numerically 

m — co ... \ 

in tardea Ila and lit), respectively. 1 Tahlesla and lb give 
t.hpir ' gr’aph i‘cal> representations. 

I % T ; * ’ ‘ * 

Remark : If the origin is moved and \ is’ replaced hy \*+<x; 

that is, setting 

£*•’.= t"+ a, l*! = £ + a, a real, 
equation (70) assumes the, form' 


1 Mr . E. Ostrow assisted with the ' computat ion of these 
t ahl es . 
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i. (*) * + JJ a rAt + AS.') = 0 

H* vest* *'t* J 


(115) 


The general f ornula ’(see (85) and (77))may now he used for 
the solutions of (115). 

Since ’ J- = -(Nf + N 3 )(sfee (7l)) if" is -found' that T a 
= J(£*-a, t*- a). Using (109) gives 

(2\*-2a) , H a = H(2\*- 2a)'. 

Thus, the generalized formula (112) becomes 

v|/(v t 9) = H [2A*(v) - 2a 3 In £g(t*> 

» 

+ V (2A*(r) _ 2a)g ( ^ (£*) (116) 

L> 2®n(n)! -! 


n = i 


£ 




j ( n+1 )(£*') s'p g^ (£*)!£*■, £* ’ = + i& 

13, The Behavior of , a "Subs on i c Flow at Infinity 


At the point -a 4 a + i3, a, 8 real, of the hodograph 
plane which corresponds to the point 1 z = <x> of the physi- 
cal plane, the stream function ^ has. a singularity. This 
fact leads to the study of , the singularities of functions 
■satisfying (46), 


If point a is a branch point, then the use of formulas 
(85) and (69)' yields -a singularity which' possesses the de- 
sired features. Indeed, 



= H ( v ) 



■£) 


T7z 


- Q. 


(i) 


( 


(a- £) 


i 



+ f « l,) 


((a- £) s/2 — a 1 / 2 £■ 


£ = A(v) — i0 


1/2 f -f a 3/fe ) + . . . i J 

(117) 


x Thl s Pi'ep ns ,!that theTyelocity ,ae. obtains 'at point 
z = 00 ' (of the physical plane). 
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is a stream function which is twp— valued at point a and 
"becomes infinite as l/(a — £ ) 1 — that is, it "behaves 

like l/(a— log v + 3,0 ) , 2 for th.e. case of an incompress- 
ible fluid . : . 

If, however, point 'a is n<5t‘ a branch point — that •• 
is, if (85) and (69) are applied to the function 
.g = l/(a-.— .£) — then 

; ■ ‘ 

X C-a - ^)’ 6 J = a Z ~ ^(y ~j + 2 ^ ^(v)(^lbg( a— £ ) ~ log a) 

+ | Q (S) (^)((c- 

is obtained, which is not a s ingle— valued function. 

For the sake of brevity the case of g = log(a — t) 

.is not ,dis, cussed, but here alsQ," in general, a many— valued 
function Is obtained. 

The function ( 118 ) can, however, be made one-valued 
by replacing the many— valued term (9 = arg t ) by its mean 
value (in the sheet under consideration). 

Clearly, this new function will no longer be an exact 
solution of equation ("46), but in many instances it will 
not differ very much from an exact solution. 1 Plainly 
this procedure may be refined. 

It is, however" ,• '•’of.' .interest from a theoretical point 
of view to determine (exact) solutions of' 5 (86) which are 
s ingl e— valued and have a logarithmic singularity at point a. 

Clearly, it is sufficient to find functions for equa- 
tion (86) which possess a logarithmic singularity. 

A function 

W»u, 8 ; X 0 , 8„) = ¥( t, li £ 0 , 

= A ( £, Is i' 0 , £ 0 > log | £ - to |+ B( i, li £ 0 . I 0 ) (119) 

1 It is noted that in this. 'case expressions (136) will 
no longer,: be complete differentials, which fact may cause 
some difficulty if a ■ is-- an infer "i : or point of the domain. 

; \ 

2 It'is assumed here-' that” 3?--’ is, replaced by 3? m . (See 

s <=>c . 9 , ) 


t )log( a— £ ) — ( a— £ )log(tx+ ! 


. . •’ I 


( 118 ) 
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which, cons‘ 14 er e cl as a function ’of , £, £, -sat isf.-i-es equation 
(86) in the whole plane, except at the point £, : = £ 0 , is 
termed a "fundamental solution" of equation (.86) with the 
af f ix at £ = £ 0 , 

Clearly, ^ = W( £, £ » a, a) represents a desired 
stream function with a singularity at £ = a. 

Notation : If 'l' — Wc , c an arbitrary constant, is regular at 

point a, the" Corresponding flow can be said to have a pseudo- 
vortex at infinity. 


■ The funct-ians A and 33 may be obtained in the follow- 
ing manner (see sec. 7 of reference 2): 


Let a new variable be introduced: 

£ = £ - £ 

... ■< * i ** * o 

Equation (86) then becomes 

d 2 ^ 


( 120 ) 




Y~ + ( £i + £o* £i + £ 0 > = 0 


( 121 ) 


A fundamental solution of. (121) with the affix at £ 1 

= 0 will be a fundamental solution of (86) with the affix 
at £ = £<°). ; 

Substituting V - a log £ x + J- a log £ 1 + B , a=*(£ ,£), 

B = B( £ , , £ 1 ) into (I2l)- gives . 

>(*£, £, + r " a ) lce K + i(~jf * Tr) 

i it*. • 

T-, +■ E ffi B = 0 * . .(122) 

Therefore a is a solution of ( 121), which has the. property 
that ^a^/£ x + a ^ ^ ;i,s r'egUlar',' at ' £ 1; > 0 ,. = 0, 
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till 


£i £] 


~ 1 ~f f ^ d ^ 2d ^ + /J 

0 O 0 0 V 


3n 


Li 


8 


/ / rm4t 3 ii 3 


Q O 


dtgdT s + ... (123) 


is a des ir ed s'odut.lqh "of' '(T2T(T.. "Aind, ’ B' i-s a ■? solu- 
tion of the equatioh' 


B 


^ 1 /+ *.» v *. o . v , *■*'■ (*?//£) « (>£,0 < 124> 


It fol-lows that _ 

,r. .. • r, - - . ■ >: £i £i * 

s - / / 0i M 1, - f f i, 


■ 0 


0 0 


f £»£ 


8 s 2 


././ 
0 0 


&d£ 3 d£ 3 


d£ s d£ s + • • • (125) 


As indicated elsewhere,. the theory of operators 
yields an alternative expression .'for a. -■ 

In references 3 and 13 a function s was considered 
r wh ich , is. w a solution of (75,) ..and , t.her efor-e which when sub- 
stituted into (76) for ,-E yields a solution of 1 (121); € 

has the following property:-. ; - , 

: ■ e (i x ,T X ,• = 1 +■ t); . (126) 

where is again a regular function of .. .£ i r £ i * (See refer- 

ence '1:2, formulas ( 1 .JL2.j , „,(i . 14) , and. (1.15).) 

If the funct ion E is denoted by e(£i»£ 1 «i ; » £o * ^ o ) 
corresponding- t o equat ion ( 131') , then 

+ 1 +1 

a < £ 1 . ) = J r e dt /( 1-t 2 ) X Y Z =? m/2-> t a I I - c 4 d.t /.( 1-t 2 ) ; 


l /s 


-l 


-1 


yields a desirepl solution,- of (l3l). .. Thus „ ■ 


a 


+ 1 

<£ , I, c„ . I 0 ) -J I-I.. 1 .’) »«/ (127) 

'• -1 ■ 


— — — — — — r--:: TT r™ — . 1 1 ' 

1 Thfe '• funct ions Q' '(v/a 0 ) which correspond to e are 

not real. This was the reason that in reference 3 and in ; the 
present paper a new solution E is introduced which yields 
-real functions Q' n '(v/a 0 ). 

In reference 12, sec. 1, the function e is determined 
In a form of an infinite series. 

It is observed that for various equations e can be 
represented in a closed form. 
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The function W satisfies •^'o'hatioc C86). Clearly 
• • dW* 


69 


= A- 


1 \ 
t 

V- 1 &o . ' - V 


. ,-vS-t 1/S 


[(X^Xo-r- V. (9 -9 0 ') 2 ;] 


+ A s .log [ ( X - X o ) 2 + (9 — 

. ■ , • ■ . ; * : : ' • 


e 0 ) 2 3 +. A. 


( 128) 


where A 1 ,A 2> and A 3 are entire functions, is also a' solu- 
tion of this equation 

w f- ! hw* 1 ■ " 

dotation. : * If • * -;0 x W*~ C 3 — — ,• C lt C £ arbitrary constants, 

o 0 

is regular at £ r a, then it can be said: that the- corre- 
sponding flow has a combined ppeudo-rortex and ' ps eudo— doublet 
at infinity. 

.... . l . . ■» .* 

• ' By re.f ini’hg thi^ ^procedure (namely, considering func- ' 
tf'ions A 4 (£< — £( 0 -)) 1 + A 4 — £( 6 )) : + " a 'b + A 5 etc.) 

other univalent solutions of ( 85 )/ bay he found which have 
singularities at .point a. . ' ' ' / 

Remark. 2; In the case where th-.e- dena itv b'r'e'ssure relation 

— 1 — , 0 iw*. 

lis of the fo'rm p - ;= A + — the functions. : K* ’ .ail'd' 

■ * 0 * 


60 


■ are 


-g- 1 0 g- 1 [. ( X — X qO ^ + ( 0 — 9 ) 2 ] . an d 


■0-9 


(X-X Q ) S + (9-0) 2 
s i l/e 


( 129) 


respectively, where X = 4- log ,i.kLZ±°l — 2 — _ — i.\ It i s 

F l[l +‘(v/a 0 ) 5 3 l/s + l J 

k 

noted that in this particular case new singularities are ob- 
tained by differentiating jW* ■ with ..r espect to X’, Bor in— 
s t ance . 


w (oi) _ 5W*_ _ ■ X - X<,, . . 

' * x ’ : ~ ’(a“- x o )' t TT0 ^T) 2 ‘ 


C 1 3 0 ) 


is a singularity which is _ inf. ip jj,.t a. of .the f irst -order and is 

aw*.- . ...■■■ .■ ' 


independent of 


■69 ■ ;• 
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14.. The Passage t.,» :the Phy§ica;l P^-ene 

Vr ~" ’ e ’ ' 

la the following a procedure will “be described for de- 
termining the-floiw in the phys deal '•pl&fc.e Qprr 6 spending to^a 
stream function, O ( v , 8 ) , giv-en ih : the hodograph plane. 


from 


Considering , as -..unknown , it is found 

50 50 50- w-i t i- r 


50 Me + Mi Mi - Mi - y Mi Me + Mi Mi = i 

bx 50 5y 50 5^. . i- bx 50 50, .. 

50 5x | 5Q 5y _ 5Q _ p Mi Mi + M. Ml - o 

5x 50 5y 50 50 ’ 5x 50 5y 50 


(131) 


;i ■ i .. 


that at' every point at which the ' Jacobian j- : 

r = , 5 C 0;o) = 50 50 _ '50 50 ' = ( y i a + v a 5 )~p = . -if -y >:• * 

5(x, y) 5x 5y 5y 5x p Q R 0 :j ' ,; ’ 

(see (24)) does not vanish and is f init e , ■ the- r elatjions 


d 5x_ 50 D 5y _ 50 

. 50 5y 1 50 5y 

D 3_x _ _ 50 *pv 5 _£ _ 50 
50 ■ 5y' 50 5x 


I; 


(133) 


hold. Using (24) yields from ( 133 ) and (132)' !; - ' '■ 

'• . ■ i*.=. dO. *. |f dO = : iO 4 r 4 ..:fL 0 . -..^: 4 i^ 4 -\d 0 

..50 50 v P ,.v 


v - U 40 + U ** * % ^ ** 

, i. \ \ : *. rY ... t! ;:-i' . : *“>, ...» u • -J 

^Since by (SO). " | " 


> (134) 


• ■ P ( 1 — M 3 ) , 

dO ="0 T dv + 'OgdQ = — — 

pv 


Oadv + ~ p^O^de • " (135) 


o. : v. 


there . is 'obtained 
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dx 


= fa/f- (1 - HfisSLi fr - g-^-1 


P U 




dy = ^ 2 - 
■ P 


V 2 "l . ' 9 

/ _ ( 1-M 5 2jBln e ^ + c os 9 ^ 


9 K- 


r in 8 


COS O^y- y Tg 


Jdv + Jj 
.■j dv + [s in — — 'l' 


1 


d9 

d6(* 


V 


(136) 


) 


or • 


dz = dx + i 


‘ « » t 

V :,r £ » 19 {.[-. * „ * i £ ] .* *[+ T * 1 % ] *.} 


■ 1 t ^(*u. e T )♦**• )(*(K«T)-i«) + T? ^ ' Mv (13T) 


Along a streamline, = constant; that is, ^ v dv + tygd9 = 0 . 

Substituting d9 = — \lf T dv /^0 into (137) gives for the 

streamlines in the physical plane the parametric representa- 
t ion*. 


x = 


y = 


r pQ c ° 3 6 [(i-M g ) % g + v^v 5 ] 

J PV E ^ Q 

p p o cos 6 l>e s + v^y 3 ] 

J pv 3 

_ j* P 0 sin 8 [(1 - M 5 )^e g + v g ^ v 3 ] 




dv 


f* P 0 cos 9 M 8 ^ I 

dv + / —2 — — dv y (138) 


pv 


pv 


T 


--dv 


-/ 


p sine Oe 3 + r ® ^ v b 3 


pv 


M/ 


dv + 


/ 


p Q s in 8 M g ^ e 


PV 


dv 


where the integration is carried out along a streamline, 

\K v, 9 ) = constant. The integrals ( 138 ) represent a general- 
ization of formulas ( 13 ) . Substituting for the expres- 

sions (69), (85), gives a parametric representation for the 
streamlines in terms of an arbitrary analytic function of 
one variable’. 


Suppose that the stream function ^ ( v , 8 ) of a com- 
pressible fluid flow is given, where 'K v , 9) is defined 

■“■It is assumed here and in the following that ty '9 and 
\J/ T are univalent functions in the domain considered. 



NACA TN No. 972 


51 




t 


4 


in the simply connected domain H (with boundary curve h) 
in the hodograph plane, and has a singularity at the point 
a (the image of z = <»). (See, for instance, fig. 2b.) 


Just as in -th : e -case of - an i-ncompr es.-aible- fluid the con- 
ditions must he determined in order that the obstacle in the 
physical plane (which is formed by the image of the boundary 
curve h of 1 H) be a c losed curve, ' Clearly, !: the necessary 
and sufficient condition in order that the image of E in 
the- physical plane be s ingle— valued is that 




+ i 



dv + 




( 129) 


wheTe l is any simple closed curve lying entirely in H + h. 


Since the integrand of (139) is a complete differential, 
t-he- value ■ of- the integral does not change if 'l ' i£ ,'con't.inu.— 
ous ly def orm ed without leaving. S + h- and without passing . 
through the singular point a. 


for 
d0 = 


Thus, in particular 
l it follows that 

- (■&. 

. \4f? 


^ &8 along 


, if the boundary curve h is chosen 

d'i' >\b 

d.\U = — dv + — d0 = 0, that is , 
dv S0 

h, since h is a streamline. If 


t h i s ._ a x pr g.s..s. laa— la -iaa-aiA e a , . . lat-a. ( 139 ) , Ih e_£a.s& i t. I p ob- 
ta ined t hat t he imag e of h in the physical plane be a , 

closed curve in a form analogous to (15) and (16). 


On the other hand, the relation ( 139 ) can also be’written 
in a different form which is often more suitable for appli— 
cat ions . , : . . 

In section 13 some standard types of s ingl e— valued 
singularities were introduced, that is, for every point a 
=‘d-+ ig functions were defined which are s ingl e— valued 
in the whole subsonic region and which satisfy eauation (46) 
there, except at point a ’where they become infinite. Such 
funct ions . are , _ . 

/ • ';W^)(v, 0; a, 0) = W*( A ,;0 ; A 0 , 0 O ) * ' '■ (140) 

. 1 1'tV is as sumed here that the ppipt " .y, = 0 . is not an in- 
terior point of 'domain H. If v= 0, 6 = 0. is an r interior 

point it is well to proceed similarly but, use' as variables v 1 
and v 2 instead of v and 0 . ‘ 
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and 


n )( v , 8; a, 3) = 


_ 5 n W ( o) (\ , 8; \ 0 , 8 a ) 


B6 n 


( 141) 


where W* Is given by (119), 3 - e 0 and 


Xo ■ j !.«r vr (142) 

4 ' L 1 + ( 1 - M 0 8 ) l/8 \1 - h( 1 - M 0 2 ) 1 A / J 


l/h. 


M 


o 



a 

-Mk - 1) a F 
2 


1/2 


N otation : With every singularity W^ n ^ , 

may he associated >a (coaplex) number X^. + 
denoted by rTw^ 11 ^. 

R(W (n >) = X n + iY n 

c 


n = 0 j If 2 1 • • « 

iY n which will "be 


( 143) 



where c is an arbitrary simple closed curve around a, the 
sense of integration being such that a always lies to the 
left. 1 


Let it be 
be represented 
in the form 


assumed now that the stream function 
in the neighborhood of a (the image 


+■ I 

r»= 7 


A n W^ n ^ + f 


i< can 

Of Z = O' ) 


(144) 


where A n are constant s and ; .if is a function which is 

regular at point a. Since curve c may be chosen for l' t 
and since the integral (13?) in which W 1 ' n ' is replaced by 

1 Since the integrand of : (143) is a complete differential, 
the value of the integral is independent of the choice of the 
path of c integration. 
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^ vanishes, it is concluded that an alternate form of the 
condition that the image of h in the -physical -plane , for a 
flow the stream function of which is given hy (144), he ~'a' 
closed curve is that 



(145) 


It will he of interest to discuss in more detail the 

evaluation of the quantities' B. Vw^ 11 M in, the special case, 
when k = — 1, ^ 


As was indicated in section 13 in this case 


W (0> = ^ log [(X-A 0 ) 2 + (3-p) s ], 


£ -P 

(x-x 0 ) 2 + (e-p - )® 


where a + i£ = a. 


X = vf- log | 


fCL-Kv/a 0 ) E ] 


i/e 


[1+ (v/a 0 ) 2 ']i /s 


:±\ x -i, 

ui ' 0 8 ° S l[l t(»/a.>j ,/a n. 


In addition to these singularities are obtained' •( in this 
particular case) hy differentiating with respect to X, the 
s ingular it ies ^ 


,(oi) = aw-C 


X - X, 


. ( X — X 0 ) 3 + ( 8 — p) 


, etc. 


( 147) 


Subs t itut ing 

■ (d)", '' 




X - X, 


[(x - X 0 ) 2 + ( 0 — P ) 2 ] v [ 1 + (v/ao) s 3 J/2 


W, 


.( o) 


.0 


■ ; C( X — X 0 ) 2 + (8 -’p )* ] 


( 146) 
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or 


Wv 


(i)__ 


: 2 (x - x o )(0 - p) ■ 


Si l/2 


[( X-X 0 ) a +(8-P) e rv[l + (v/a 0 )*] 

w (i) ( X — X 0 ) a - (9-- fr) s 

8 " C(x- x 0 ) B + (e - p) 2 ] E 


(148) 


or 


w ( 01 ) = 


— (X — X(j)? + (8 — 0) 


[<x-x 0 ) s + (0-p) 2 ] 8 V [l + (v/a D ) 2 ] 


2 ll/2 


y ( o I ) , ~2'(^"X 0 )(9-p) 

’ e C(x - x 0 ) 8 + ( 6 - p) 2 3 s 

into (143) gives the corresponding values R^W^)^ = Xjj + iYjr, 
k = 0,1, and R^W^ 01 ^ = X 01 + iY 01 , respectively. 


For instance, 

X 0 = /. 


jT -COS' 

IL V 2(1 + 


0(9 0 ) 


s in 0( X — X 0 ) 


2 \i/ 2 


[( x -x 0 r + (8-e) ] 

• + j^ CQS 9 (X-Xq) + ( 1 4- ( y/a 0 ) " ) 3 ( 9 - P ) s in e j 


dv 



p{E 


— sin 6( 9 p ) 


cos 0(X— Xo) 


[( X-X 0 ) 3 .+ ( 0-0) 3 )- l- U a ( 1 + Cv/a 0 )*) U 




dv 


+ 




sin 9 ( X— q ) 
v 


(1 + ( v/a Q ) g ) 1 /s (8 — 0) cos 8 
v 



■ 15-.- Appendix -■ 

I. Verifi ca tion that th e expr essions ( 156 ) are compl ete 

d iffe re nt ia- J.8-.-~ To show that dx as given hyi (136a) is a 

2 By ' (lS6a ) -'wni 'i>e denoted in this section the first 

expression of (156), and by (136b) the second expression. 
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complete differential, it • ds ' iieces s’ary to'prove that.' the 
coefficients a and "b 


a = P, 


[- 


p.Q. g . gii-- 


= P e _ 


py 2 
cos 9 


- 


s in 6 
pv 


* 


>3 




(149) 


of dy and d0 of" the right-hand side of (l36a) satisfy 
the relation 


5a _ 5ja 
58 5y 


( 150) 


Recalling (30.) yields 


5<t> _ pQ 7 ^ 5<t> p 0 ( 1 - M ) 5\j/ 

50 p 5y ’ 5v pv 56 


(151) 


s o. that 


cos 8 - ^ s' i'n 

TT ' ■ -v 


Ther ef or e 

5a _ 

‘5V “ p o 


• a = <=o ( 

» - Po !£*$*■*»' -•£**») 


£-I'*S ' 

pv 7 / • 


( 


sin 8 ^ + eos'9 ' ^ . _ cos 8 


V0 


\k - 
v y 


s 

s id 8 * 
py 


* 


(152) 

ye) < 153 ) 


5b / « c 

■5? ” p 0 ^ 


cos 8 

— 


<jy Q + cos 9 ' * irf 8 p+ 


0y 


/ sin 8 ,,, 

+ — Vq - 


py 


P 

sin 8 
py ■ 


II \lr 

2 y v 8 


* 


T ®) ‘ 


(154.) 


Thus it ’is necessary to 'proye that 

/sln^. , _ = _ M. e’ . - ’ _ i m. . ... 

V • T . T. . p- T V .... . T » .6 p*. . T. 

. .. ; '- ; + «iajS.' 




- .• V 

> - . - 


(155) 


• py 
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If use is made of (30) again, it I s ' neeeb sary only to show 
that •* 1 < 


I s: 


• • 1 - M‘ 


pv’ 

,2 


p -v, pv' 


' M = ,-v ■ (log p) 

..But from (25) it is known that 


P = P, 


k - I v 3 

2a 0 s 


k-i 


so that 


v -iL (log p) = - 


SV ' ' a " - I (k - .l)v* 


(156) 

(157) 


(158) 


If this is compared with (3l) it is seen that this is exactly 
equal to — M® and it therefore has been verified that the 
right-hand side of (l36a)„, is a complete differential. In a 
similar fashion it might be shown that the right-hand side 
of (136b) is also a complete differential. Hence, since this 
is the case', it follows that the integrals (138) are inde- 
pendent of the path of integration. 


II. A. .proof of an auxiliary lemma .'— In the following 
it will be proved that F(2X)(see sec. 9) can be approxi- 

m 

mated by polynomials F^( 2X ) = ^ m ^ e® s ^ , a*' 12 ^ = 0 

/ ■ s o 

S = 0 

in every interval (- 00 -, X 0 ),' X 0 < 0, and indicated how to 
determine the F n (2X). 


a 

X D < 1, 


If 2X is replaced by log- 1 X,-: the F[2X(X)] is 

continuous function of X in the interval (0, X 0 ), 

and by classical results , it. is obvious that it can be ap- 
proximated by a polynomial in X. It will be' seen t.hat it 
will not suffice merely to approximate F, but in addition 
to this it will be possible to require that any given number 
of the derivatives of F be appr oxima.t ed in the interval 
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It is, however, of intar.egt tp give a more explicit 
form of th : e a§>pT'bxia§,t ing.. polyno’mi.'als . • tPh. is w ilil.‘ make it 
possible t o \det ef.mine th‘e*. corf.ect ionajof, obtained 

in section 1? which have to be made in order to obtain func- 
tions Q,( 11 / corresponding to a given P m . 


i *• | J,. , >\ ■ 




: r. ! L 


If M increases steadily from 0 to 1, X increases 
from —o° to 0. Since the relation M — s» X is a one— to— one 
correspondence, t^ere cor.r,esppnds to every X 0 , X 0 < 0, an 
M 0 = M( X Q ) < 1 . Therefore, if X e( — ® , \ 0 ), then M(X)e(0,M o ). 
F or M. = M 0 , the funct ion 


j. _ (k 4- 1) M 4 [~(3k - l) M 4 - 4(3 - 2k) M 3 + 16] 

64(1 - M 3 ) 3 

may be approximated by a polynomial F n . of the ( 2n + 8) 
degree in M 2 , 

F n (M) = --~ A - 1 ' ) M 4 J^ — (3k - 1) M 4 - 4(3 - 2k) M 2 


th ■ 


64 


n 

+ is] r y (b- 3 )(-i) u h s " 


i 


J L 


U=1 


(159) 


Only a 'finite number of powers of M 3 appears, in ( 149 ) . . It 
will now be shown that ■ .M 3 can be developed in the uniformly 
convergent series- • : 

U 2 _ V a v K v _ „ f ( k+1 ) 1 /3 - ( k-1 ) ^*VE=x s x 

M _ / p K X * x ", 5 f: T7S — • r~jr) * * e . X < 0 _ ( 160) 

^ ■ - A(. k+ i) 1 /a + (k _i)*/V . , . - 


Instead of considering-* M?-, it is well, to introduce 

s = 1-(1-M 2 ) 1/2 


( 161 ) 


Since M 3 - 1 — ( 1 — s ) 2 , it will suffice to determine the 
s er ies for s . 


From (48) follows 
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X 


S . 
2 — S 


^ (•h^ 1 -+l'+ S ; )- ' V.-.L+ fc- l .< ... > ■ " 


■ • h'- 



k > 1 


(162) 


For simplicity’s sake it will be assumed' in the following 
that 





k ^ 5-/ 3 


(163) 


Bow consider the function X = X( s ) . as a function of the 
com p lex variable s, and investigate its behavior in the 
domain [ s j < 1 . 


Remark : X( s ) = — ( — — ) h ^ is a many— 

2 - 8 V 1 - 1 + b / ' ' 

valued function because any integer may be taken for H. 

Since, however, its branch points s^l + H -1 and -s = h“ l — 1 are 
outside I s.J <1 it is necessary only' to consider p.’he of its 
branches; Therefore, H »• 6 ""is "chosen, so that whenever 
X( s ) is mentioned this branch will be always understood. 

\ { v r . * 

In order to pr.oye 'that tjhe ‘image of isj^l is a schlicht 
domain in the X— plane it is' noted at first that X( s ) , S 
real, is a real function,, and .therefore the image will be a 
domain which is symmetric with respec-t to the real axis. The 
im.ag^ of,, s =• + 1 will be' the point 1, and the image of s 
= —1 will be a point of the negative real axis. 

It will be ghown that,- If- •'<!> = a'rg s, varies from 0 
to rr, X-increases s t e a d il y~. ' ET^trti n g s = yields 


1X1= ~T _ - _ 

(5-4 COs^) l/2 


(l + h- 1 ) S + 1-2(1 +h Vost? (164) 
( h -1 — l)® + 1 + 2 ( h —1 — 1 ) cos cp 


i /2h 
M 
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and 


d[ X 
dct>~ 


/ 4 sin 4) \ j .1 

V( 5-4 = o s 4>) 3/ = ) 


/sh 


! 


+ ri/" — i — 

L \ ( 5—4 coa <P) l/3 J -1 

f (l— h~ 1 ) S +l+2(lT 1 — l) cos <£>")( l+h~ sin<}>+3(h~ 1 -l)33ln0( (l+h~ 1 ) e +l-3(l+h~ ^~) C OS 4?) 

^(l-h _l ) s +l + 2(h -1 -l)cosc^ S 

= ( ' s - si ^ - /a ) - : : — P (165) 

V(5-4 cos $) ' (l-h ) +l+2(h -])cos $ 

p = ; f- _ r ( l+h~ 1 ) 3 +l-2( l+h~ *) cos 4> 1 - • 

L \ 5 — 4 cos 0 j 1 ’!■ 


— 4 


(h _1 -l) 2 +l+2(h 1 -i)cos $■ 


■J 


= 4(h 3 -l) 


(h +l),— (h s +2)cos 4* + cos 5 0 
( 5—4 cos 0)^(h 1 -l)‘ 2 +l+2(h ^ljcos 


For:, all values' O' ^ $ = • n the: expr??§ion (165) ha.s the same 
sign as P. The denominator of P is always positive, and 
the numerator is positive for — 1< cos4><1. Clearly for all 
values of 4> , 0 < 0< tt, P , and .therefore ( 16.5 ) is _ pos itlve . 

Thus the "boundary curve of the image of | s } < 1 .is a~. 
curve which does not intersect itself. By classical theorems 
of the?, theory 'of functions the domain bounded -by. this curve 
is schlicht; Clearly it includes in its interior the domain 

[ Xt < 1 . 


Since the image of laj < 1. is schlicht and includes 
IH<1. the -inverse function s = a(rX) is regular in 
1*1<1 and by Cauchy’s theorem can be expanded in [X} <1 


in the f‘or,m of an infinite series 


s(X) 


- -V, 




v 


v = 1 


For" every X D < 1 and every e ’> 0 ■ there exists an F 
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N 


N 


such that j s(x) — V"* Pi>X U j^£ for [X|<X 0 . Thus ^ Pu^ v 

• U = ! . I 


y i elds /.-th e required approximation. 

.t 

Remark : Clearly- N can he determined so large that any' 

N 

given number of derivatives of ^ j3 v X U approximates the 

corresponding derivatives of s(X). 

It is noted further that a formal computation yields 
for the right-hand side of (152) for [Xj < 1 


X = .3 


- 4" (3k + l)s* + ~ (4k® + 2k-l) 


4 


- ~( 24k 3 + : fefe s -l4k-l)s 4 + 48k 4 +4k 3 ~44k a + 2k+5 ) s £ 

,24 48 


48 0k 5 -104k 4 -572k 3 + 14S.k s +l 2.6k-25 )s® + -g£~( 2880k® i 

-1584k 5 -3944k 4 +2212k 3 + 1140k s - S02k - 5 ) s 7 + . . . (166) 

The inverse function is ,i. . - 

b '=-'X’+ -|- ( 2k + l)X® '+<. .^415* + 6k + 3)X 3 +'i-(24k 3 

. ; „ J - + 68k s U76k + 29>X 4 + — (48k 4 +212k 3 + 392k S + 328k’ 

•+ 103 )X 5 + — i— ( 480k 5 . + ' 2976k 4 + 7968k 3 

48 0 ■ • ■ • . : 

+ 10 7 8 8 k S + 7266k + 1935 )X® + (2880k* 

+ 234.72k 5 + 84232k 4 + i'62124'k 3 '+’ l?394fe'k'® ’ : 

+ 98086k + 22675 ) X 7 + , . . (167) 

By the present result "this series converges-, in- [X]<1, 
therefore for 0 c X < 1. 


and 
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•• 1 III. ADDITIONAL REMARKS ' 

.1 

16:. - The Boundary Value Problem in the Physical Plane. 

Mixed Plow 

• i ' 

The theory developed in the second part of the paper 
leads to various methods for constructing flows around air- 
■£ ci jL.s . 

The primary problem to be faced in the theory of airfoils 
is to determine the flow with a certain velocity at infinity 
around an obstacle given in the physical plane. This leads to 
a very complicated nonlinear problem in the hodograph plane 
since the domain where the flow is defined is determined by 
the flow itself. However, this problem may be considerably 
simplified if it is agreed to obtain a flow around an obstacle 
which approximates the given obstacle. 

The hodograuhs of flows of an incompressible fluid around 
prdfiles of certain types and for a number of angles of attack 
may be determined once and for all. 

The present approach also makes it possible to construct 
functions satisfying (3-2) and having singularities of the 
kind required - that is., singularities of the flow of a com- 
pressible fluid whibh yield sources, vortices, and doublets. 

A hodograph is chosen which in the case of an incompress- 
ible fluid leads to the desired profile. 

Let $(v,9) be some solution of (32) which possesses 
the required singularity at point a (the image of z = oo). 

A solution of \J/(v,6) of (32) is further determined, which 
is regular in the domain. H and such that 

. $( v,G ) + \{/( v,6 ) •• 

assu-m-es a constant value on -the boundary h of H. 

'-The obtained function is- ; a hodograph of a flow of a com- 
pressible fliiid 1 the ima'ge o-f Whi’-oh in the physical plane will 
in many instances not differ considerably from the given pro- 
file. Thi g 'method' -of a : t : tack 'can be refined:' By the forego- 
ing procedure the initial profile is distorted ih a’ certain 
way; if the given profile is distorted in opposite directions 
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and if the procedure described is repeated to the distorted 
profile, then in many in st an.qe s ; a; fbet t er ; #ppr oximat i on is ob- 
tained. This method may be repeated un‘ti’1 the desired degree 
of accuracy i s at tained . 1 * * * S 

However, this procedure h'as the inconvenience that in 
order to determine \|/ it is necessary (at each step) to 
solve a boundary value problem for the equation (32) which 
requires rather long computati ojr. :In another paper the 
author has developed' in ‘de-bai-l- -an ■ alternative t o .this ‘ method ,, 
in which he avoids the necessity of solving boundary value ‘ 
problems. . . f 

In the present considerations' attention was in the main 
directed toward the subsonic case. In addition to the method 
of attack, which is based on considerations of section 8 of r 
the second patt , ther e exi st s another possibility for handling 
the mixed problem - that is, to construct flows which are'-paf- 
tlally subsonic and partially supersonic. . 

17. The Representation of the Stream Function of a 

Subsonic Flow in the Region in Which .the Velocity 
is Near the- .Velocity of Sound . 

V ...» 

Part ial ly -.Super soni c Flow 
* - , 

In the region < M Q < l) where M 0 is near 1, the 

series (85) converges very sldwly, and it is .therefore neces- 
sary to employ a large number of 'terms in order to obtain a. 
good approximation for '4/*. If this be the case,, it is then 
expedient to replace the expansion (85) by (103)". 

1 ’ ■ - l .... . 

This is, however’^ not the only' : way of 'overcoming this . 
difficulty, and in the following, other means of.' so doing 
will be indicated; this alternate approach employs the method 
of "analytic c ont inuat i on ; , , 


1 lt may be observed that a similar procedure can be ap- 

plied to prove that for every profile (satisfying certain 

conditions) there exists a flow of a compressible fluid. 

' . ' r. t 

S " 

This method will. -be. developed in more detail ,in a fu- 
ture reoort of the ah.th.or..-. 
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Let ty(v,Q) 
|\f/ n J,n = 1 , 2, . 

iut£jQ.H.:S of (46), 

t * 

Suppo'se that H 


* ' . *r *' 4 ' _ 

’be detefmintfd/'itt a domain, say H, and let 
, be a "complete" system of particular so- 

each \}/ n being det er mined ■ in a domain G. 
and -G- actually do overlap and denote their 

03 


common part 
pansion of 


by I. Further, let 


y a n 




n 


be the series ei- 




\J/ in I. Frequently ^ a n^n v wil1 converge 


n=x 


out side . of 


I, ■ ^ay in the domain H a - I- , where H 2 is G 


or some par't of it. If, in addition, 

wise differentiated twice in H 3 , it 
continuation of \|/ in H 3 - I. 


00 

I can b e t er in- 
ns'! 

represents the analytic 


Remark : The requirement that ) a h'i ; n coincide with ' in 

n^i 

a domain I, can be replaced by another requirement, which 
will be explained later. 

Freauently, the dpmaip. H a in which the stream function 

* ' • . v Co * s » ... ■• . 4 

can be represented in the form y s-n^n covers a supersonic 

r n=i 

% 

region as well, and- cphs-equently. thi s method will then. yield 
the flow in this latter region.- In ‘this manner, a method 
(based on considerations other than those of sec. 8) for de- 
termining a mixed flow may be obtained. 

Two alternate forms of this method will be ‘discussed in 
the following. 

* ‘ ’ • * . * t 
First Method 


In order to develop the first approach, an auxiliary 
lemma must first be proved. 


Lemma : 
an analyt i c 


Let p(v,6), v 0 < v < , 

funct i on of two L real- variables 


L < 6 < L ’ be 
v , 0 , and ret 
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CO 

^ a 1) (v) cos 2^- :: ^-'' , 'b u (v). sin 


-rru6 


u=o 


a 0 (v) = — L*" 1 p(v,6) d0, a' (v) = — P p(v,6). co.s. — d0 , 

2 / L / L 


-L , 


-1 


Id (v) = — / p(v,£.). sin d0 

^ L / L 

-L 


(168) 


be it s Fourier development . The series (160) converge s unl - 
f or mlv and can he di f f Sr ent iat ed t er mvi sa any finite number 
of times both with r e snsct t o v ' and with r e sue ct t o 0 . 


Proof; Let 


S^Cv.iQ ) 


. 3v 


k 


Z r, k 

d v 

• u 


d ( v ) Tfu0 d K by(v) 


cos 
k L 


dv 


k 


sin 


TTU0 


.Now., since p is an analytic function of v and 0, 
is also an analytic function, and therefore 


, ( k e 1,2) 

■S 3 P 


dv 00 



a 3 p 

. ^0 V 00 




dS < A 


(169) 


-L 


is bounded, uniformly, in v; . 
On the other hand,. 
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from which the uniform convergence of the series 

CO 

T 


da (v) tto 6 db ( v) 

— 12 cos + - — U sin 


(171) 


dv 


dv 


follows. • But (l7-l) is the series which, is obtained by -'differ- 
entiating (168) term by term. -In a similar way the other j,... 
cases may be handled. 


Since eyery - solution of an elliptic equatign ,with ana-.. 
lytic' coefficients'-i's an analytic idaetl on^of-two real-Vayi— ■ 
ables, the result obtained cai be apnliedto the case where 
p .( v , 6 ) is the stream function \j/(v,Q) of a subsonic flow. 
Thus 


' ’do ’ 

.. X 


U = 0 


a ( v ) cos H2JL + b ( v ) sin - 

u ” Jj 


a 


jj 

(v) = — / i f(v,e) de , 

2L -X 


ay ( v) = J- / \i/(y,0 ) CO 

•J.L 


s Hi d e, 


* v-.-. ! 


■, t-’: f 


TV- 


d p U) = y \i/(y , e ) sin HHi de, (u = 1,2, . . .) 


-L 


J 


can be differentiated termwise. If now, following Chaplygin 
the author introduces Instead of v the variable 
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T = 


v s . ; : 

2a o' 3 


then the ea.uation for, , \j/ assumes the form 


^[ 21(1 - + ± - r-J *L - ± .11 (1 _ T)“ g = 0 (173) 


BrJ 


2t(1 ~ t) 


ae a 


where p = 


(k - 1) 

Differentiating termwise gives 


(reference 1, p. 5, formula ( 1 2 ) ) . 


y (A{t<i . vv- r i£ky_ j - ye ± y t u . T v- e aS^ ieo. 

L IcLtV dT J T.( 1 — T ) ■ 4.L- 2 L 

u=o L 


(t(1 - T )~ P - Alii (1 - t)“ P ulnf. h i sin 

V ,J dT J- . , T ( 1 - T ) 4L 2 V J, 


d_ .4 

dT 


utt6 

L 

(174 ) 


and, therefore, the a., and h are each solutions of equa- 
tion ‘ V ' 


A 

dT 


A 

dT 


fa - T)" e ^4 . 1 - <38 + 1> T (1 _ T)- s ' IJLS . 0 ) 

l dT / T (1 - T ) 4L 2 


> (175) 


(V ( 1 - T)" e iV) - 1 “ (2P + 111 (1 _ T)" P ” S v\ = 0 

V AT / T (1 - T ) : 4L 2 


(175) is a hype rge ometr i c series and thus e^cry solution of 
(175) may he written in the form 


JL. 

. . 3L 

.l*v*v + , (176 ^ 
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where » and !?B V . ar'fc const'ant-S and 

= Ftoy.py? -6; 1 - T>, 

; "" • i ’* ' • • ‘ ■■ r 

V = (1 - T)" p+1 F( Yy - ay, Yy - Py; 2 + P; 1 - 

* r * " . i • . 

„ ■ ' V f -i 

P(a,P,Y,T) being the hyper ge ometrj.c series. Here 



In order to determine the constants Ay, By, the 
theorem is employed: . . _ 

; - » 

- ( l ) . . , ( a ) . / . I : I 

Let v|/ - ... ( v,8 ) and \|f . .('▼> 0 )l [ y 0 < v < v x , f-L 

“be solutions of an eauat i on of elll-ptic type : ' If '. 

line , say (v = v 0 ) 


’..TP 


* 


(i) -(y t (s) ( v , s) 


and 


a^ (l) (y,6) 

By 


y=v 0 


9t (3) (y,8) 
By 


Y sY o 


then i n the whol e domain [ y q < v< , - L <. .Q.. < 


■^■ l) ( y, 6 ) = (v. t ,e ) . 


\ (177) 
T) [ 

- 

following 

< G < L] 

' "along a. 

(179) 

(180) 

L ] 


(181) 
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l. . ■ 

Remark : Suppose that the function g(£), £ ~ X ~ i8i 

X = X( v ) (eee (85)) is regular in some domain Hj +. H s , 


H i ~ |^ v o < v < Yi» - L < 8 < L 
H g = £v x < v < v a , - L < 6 < L 

3 3 

which domain lies in [6 < 3'X f X < 0], Then by the main 

theorem it follows that ^ = I,m.P(g) is also regular in 
Hi + H s . Suppose that \|/ has been evaluated in the domain 
Hi, but it is desired to avoid the evaluation of by 

means of (8.5) .■since this series converge &- very slowly in H z . 




% 


\ 


CO 

y S (l) (t) cob 221 + s (S > ( T j) sin 221 (182) 

■ v=l 

UL 

( k ) f ( k ) • ( k ) , -] 2 

'S v (T) = U v + By T , (k = 1,2) (183) 
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he'the solution- \1/’ under consideration, the constants 

(k) (k) 

, B p must he determined so that ■ 


<s ) f \ 

s u (.T >. 


S„ (4) (t> 




V=V, 


V=T 0 


= S ( ^0 ) 


= 


V 


, i- 


- (T) 

dv 


v=v. 


d a ^ ( v ) 
dv 


v=v, 


d8.„ (a) (t) 

dv 


v=v f 


dh^v) 

dv 


v=v 0 


It is notioed that 


dS 


(k) 


dT 


A (k) (k) 

u dr u dr 


. J2_’ 
• ■ 2L 


v 

V 2L 
+ — T 
2L 


-l 


A 


. (k) y + b P 


V V u 


-•] 


aF ^ = + 1, e*, + l; -P + l; l - T) 


d ( 1 ' — : T ) -S 


dr. * 


d ( 1 - t) 


-•= (1 + p')(l - T) p ; p('V uT a v ;.Vy - 6 V ; 2 + |3; ;1 


+• d - t) p+1 ~ ..gttj ( , Z , p ,P_ p(v ■-■o u ‘+ iv 

2 + P 


Y v - P v + i; 3+f3; 1-T 


(184) 


(185) 


(186) 


- T) 

■ , M P 


4 
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Since in 


v ; < v < v 3 , -L < 6 < L 


function, the series 
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C t - , • » 

\Kv, 8) is an analytic 


f {[V 0 T « + B » <0 V 

.u=l 


uL 

2 


cos 


. ttp6 


.uL. 


A ^ s ^ F + B F »1 2 sin HHif l 

u uu UJ 


(187) 


ani its derivatives converge uniformly and absolutely in this 
domain. Then (187) represents the solution \|/ under consid- 
eration in the region H a . Moreover, this series (and its 
derivatives) may also converge outside of R z , say in 

H 3 = [v g < v < v 3 , -L < 0 < L]. If H 3 partially lies out- 
side of the domain [9 s < 37\ s , \ < 0] (see sec. ll) then the 

obtained expression gives the analytic continuation of the 
solution outside of the domain of representation by the inte- 
gral formulas (85). In particular, H 3 may include some 
region which lies in M > 1. 


Very often it is known that the region, say L, where 
the velocity is supersonic is small. Now, instead' of summing 
t o inf inity , take 


N 

y [s u (l) ( T ) cos + S v (s) (T) sin J2jfi 

v=i 


(188) 


(see (183) and (177)) where N is sufficiently large; then 
(188) can be considered a sufficiently good approximation for 
analytic continuation of the stream function \{/ . under con- 
sideration, On the othVr 'hand, (l>88) : r e?pr-e sen-t s ^ in. the 
whole plane and therefore is particular in L. 


In this way are. obtained .approximate flow patterns .which 
are partially super sonic. 1 Iii applying this me.thod,. it’ is 
necessary, however, to check whether the streamlines in L 
approach to smooth limit line ^ when m increases. 
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* . .. , r . • i Sec oad'. Method - ■ ■> 

_■ . ..In reference. ,2 the autb©,r ; ,-ha s introduced ’.differ erit' - " 

■•'methods for computing sets of f part icu\ar solut i one, 1 jj/ n ( v ,0.) , 

Of (32). (See ,'p . 17 arid p. 2^ ’-of -reference .2*) i ;.T.he fu.nc- 
* tions of each Of these sets are defined for the subsonic and 

the supersonic range. 

Let H "be a domain in which it is desired to determine 
a hodograph with a : super sonic velocity. Then H is divided 
into two overlapping , parts Ej. and H g . .In,.. Ej the veloc- 
ity is throughout subsonic. The intersection of H x and 
H s is d.e.not-ed by 1 I. ■ In figure '5, Hj. 'is' that part of H 
for which v<v x , and H s is that part for which v ps v Q ^ 
v Q < v x < 1. There is determined a function \J/ 0 (6,v) which 
is defined in • Ej. an,d.-hae at point a, a prescribed singu- 
larity, and on the part of h which lies in v <'v lt h 1 , 
approximately constant values. 

r " ■ How' consider the' functions 

/ n 

‘ 5 ’ ‘ \k p (v.,G) + y a v \|/ v (v,e) 

" • ' ' '■ ' '' ; v=i 

m 

. . V*? = .y • . 

V=x ■ ■ 

and determine the ct^ and in such a way that 


n 

r- — i 

■ .• V ' p 

.* ’• n .. . 

^o + y 

‘ is + . J ., 

y Vv. 

U=1 

rt ■ 1 

■' ^3 

V=x' 


ds 


I" • • ' . ' v-i' ’ 

will be a nlinimum. 


-fi, t : . 3 . . ■ *. . 

+ y <*»% - y % dvd6 


U=1 


r , i > -7 


It is observed that it is possible also to use the 
Chaplygin solutions. See, for instance, reference 2, pp . 18- 
22 . 
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% • *4 » > * V I m r * 1 1 * 

If the boundary value problem' has a solution (possessing 

certain ’pr oper t be S' bn .the. "boundary ' h) and. if. the syst.e.m. . "f.ll/y 

i's' complete , th.en it is possible to show, (under, .certain, ad-di- 
,ti”onal conditions) that the limit .function obtained by thiis' 
‘proce'ss will yield the solut ion, ■ 


....... 1.8,, A R'emajp-k Concerning .tjhe Applicat ipn>-.of .the. • 

Holograph Method in the' Three-Dimensional Case .« 

: .The- method, developed in this paper yields a general 
formula for the -stream functions of possible compressible' 1 

fluid flow pattern^. 

\ . r- y • 

• '■ As indicated in reference 2 (secs.’ 6 to' 8) there exist 
other methods of obtaining particular solutions of equation' 
(32),‘ and for deriving from them solutions of- (117). They- ' 
often are not very conveniont for practical purposes, and in 
many instances represent a flow only in a part , of its domain 
of definition. 

In the following will be ., i ndi cat ed a method of obtaining 
particular solutions -of equation (32) which has the disadvan- 
tages indicated but which can also be applied in the three- 
dimensional case. 

As is well known, the Velocity q = (u,v) of an irrota- 
tional fluid flow satisfies the equations 

V (p if) = 0, V x q = 0 , ( Cauchy r Riemanp equation) (189) 

* ■ ■ * •'•* v , V * 1 ■ ■* ■ •» .■ 

This suggests considering three-dimensional flows where the 
velocity 5 = (-u,-*v,-w) satisfies the equation 

V(p q) = 0, r V x q = 0 (190) 

„ 1 

p = p(V) being function of V = (u 3 + vf + w 3 )^ alone. 

It follows from the second equation of (19(5) that there 
exists a potential <p , such that . 

■ r 

q = 


'(191) 
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Inserting this'Value in the first equation of (190) yields 


8(p dcfr/dx) + d ( p 80/dy) + 8(p a<$>/8z) _ 
dx By * ■ di 


(192) 


(192) is a very complicated nonlinear partial differential 
equation. , r . t ’ . 

The introduction of u, v, w as new variables leads to 
a much simpler nonlinear differential equation. 


Introduce as new variable^ 


u = 


d(p 


v = 


d<P 


w = 


8x 4 _ dy 
and as the new unknown function 

X = xu’ + yv + zv - '4> 
Use (193) to obtain from (192) 


* i_ w 

3$ 

3z .. *>•' 


(193) 


(194) 


8(pu) + 3(pv) d(pw) 


dx 


ay 


3z 


p u + p 

u ... r 


8u 


, 3v 


+ p u — — + p u 
dx v 


3w 
w” 3x 


+ p„v #£• + 


ay 


Pv 7 + p 


9y Pw By Pu 3* Pv Sz 


, + 


[rw» * p] - 0 < 1S6 > 


It follows from (194) that 


X = 
u 


gx ay 3z 

x + u2_ + vf^-+w^- 
8u 3u du 


ax ay 34> az 

dx 3u ~ dy 3u Bz 3u 


~\ 

x 


X^r — y I 


K = z 


r d96) 


and 
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ax 

X 

Bx 

% 

Bx 


Bu - ' 

uu * 

" 3v 

f ll'frt 

' 3w 

Vw 

By 

A __ _ _ . 

a z 

- X 

B z 

x 

3w 

vw • 

Bu 

UW' 

3v 

vw * 


From . . 


By , x By 

a* : ' . U Y ! .av 


a 2 * 

= A WW> 


= 




— = — x + — x + — 

au ax A uu s y A uv a z 


Niw - 


= i 


au _ au 
bv ax 


uv 


+ ^.A 


8u _ 3u ^ 


By 

3u' 


■ vv 


A. 


+ ^ X 
3z vw 

3u 


= 0 


3w dx uw 3y vw az 
there is obtained 


Vw = 0 


Bu 

Bx 


Vrw 



;D, and so forth 


) (197) 




> (198) 


J 


(199) 



■t- 

l . 

- 


the 

determinant 




x ;uu 

Vv 

^ uw 

' D 


Xv 

Vv 

Vw 



Vw 

^ vw 

^ ww 


Substituting the values obtained in (199) into (195) yields 
he following equation for A: 
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X 

X 


X 

X 


x 

X 

[pu u + P 

* • * , " r 

vv 

vw 

- P v u 

uv- 

• vw 

• 

+ p u 
p w 

A uv 

A vv 

^vw 

^ww 


^UW 

A ww 

^uw. 

X vw 


- Pu v 


.+ P u w 

t 


^UT ^UW 




^UT ^UW 



r ‘ "! . 1 

*uu 

A uv 

' «.« 

A uu 

A uv 

+ 

Pv T + P 

^wu 

A ww 

Pw v 

A ut 

A vv 


. ■ 




• 



- P T V 


‘ vv . "vw. 1 . 


^uu \v 


^uv ^ww 


+ (p w w +p) 


^uu ^uv 


^yiv ^yv 


2 2 2 2 

Her.e .p- = p(u + v + w ) is a known funation; , 

i -.There now arises the problem' of determining particular 
solutions of (200). Clearly, this can be dohe by fusing the 
series developments 


/ , man 


■ : (201) 


m , n , p . 

which satisfy equation (200). ' ’ , 

Such a series development which represents' (in the hodo- 
graph space) the potential function (£> of a possible flow 
pattern of a compressible fluid converges! oniy ip the neigh- 
borhood of "the origin. . ! ' ■ 

However, there exist method8 v of determining <t>, in the 
whole region of the real (u,v,w) space where cf> is regular. 
Such a 'representation, for instance, is given in many cases 

tty . ‘ ■ . ; • ' • 

4 >(*. 7 .«) - iim ••• y 

k ^ 0 * — r[l + k(m + n+p)3 

m,n,p 
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OONO-LTJDING- REMAHKS ■ • . ' j ■ 

’ ' ■ ‘ | .- ( 't • . ! ; . 

The main result of the present report consists in deriv- 
i ng va- ’’•f ormula which- transforms an arbitrary analytic function 
of: a complex .variable into a stream function of a compressi- 
ble sub $on'i'c flow. . 

This formula yields compressible flows around symmetric 
(and certain ponsymmet ri c ) obstacles. 

T,he' main difficulty arises in adapting the formula to a 
given shape of the obstacle. Approximate methods for solving 
this problem are indicated in section 15. 

Since all expressions appearing in the theory of a com- 
pressible fluid flow are much more complicated than those 
occurring in the study of incompressible flows, a careful 
investigation of the numerical methods to be applied is nec- 
essary. 

...... A considerable part of the numerical work consists in 

preparing tables of auxiliary functions such as Q' n ', which 

have to be used in all particular ..cases. In this paper the 
functions are computed up t o n = 4, for k = 1,4. 

Tables for the for higher values of the- super script 

n will be necessary if flows with maximum Mach number ap- 
proaching 1 are to be considered. , . 

‘ Bach particular problem also involves the- performance 
of certain integration processes. In order to advance the- 
application of this theory it would be necessary to use ef- 
ficient modern ' cCiaputihg' devices. ■ 1 .* • - 

N* . • r I , * • 

■ The present paper ' dpals only with subsonic flows. It, * 
should be emphasized that the development of the theory will 
permit consideration of flows for which the maximal velocity 
exceeds that of sound. (See sec. 1?.) -• 


Brown University, 

Providence, E. I., May 15, 1944, 
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Table It 


Values or F, H, -q<T~° 
Q 0 * 1 &"!,*] roK k-/.4 




mmmm 

Hi 
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A 70S 
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Table Ila 

The values of F, H, Q^for 


2A 

u 

T 

v/a 0 

H(2X) 

- 

0 

1 

0 

1 

-2.82 

0.265 

0.964 

0.272 

1.000 

-1.96 

0.403 

0.915 

0.430 

1.001 

-1.612 

0.473 

0.881 

0.518 

1.002 

-1.51 

0.497 

0.867 

0.551 

1.003 

-1.388 

0.528 

0.849 

0.593 

1.004 

-0.79 

0.690 

0.723 

0.861 

1.015 

-0.394 

0.820 

0.572 

1.165 

1.046 

-0.188 

0.901 

0.433 

1.440 

1.111 

-0.132 

0.921 

0.388 

1.527 

1.145 


Table OT>. Tha values of F, H, for k = 1.4 


2> 

■9 


F 

H 

Q C1> 

q (2) 

q <3) 

-«o 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

-j.ernz 

.1000 

.0999 

.0001 

1.0000 

.0000 

.0001 

-.0001 

-2.5096 

.2000 

.1992 

.0011 

1.0002 

-.0009 

.0015 

-.0013 

-1.7327 

-3000 

.2972 

.0064 

1.0014 

-.0056 

.0085 

* -.0083 

-1.2071 

.4000 

.3938 

.0256 

1.0042 

-.0199 

.0342 

-.0395 

-.8238 

' .5000 

.4879 

.0866 

1.0110 

—0574 

.1160 

-.1680 

—6706 

.5500 

-5341 

.1565 

1.0167 

-.0935 

.2101 

-.3496 

—5364 

.6000 

.5795 

.2839 

1.0247 

-.1501 

.3823 

-.7494 

-.4204 

.6500 

.6242 

.5245 

1.0359 

-.2401 

.7089 

-1.6847 

-.3203 

.7000 

.6680 

1.0060 

1.0515 

-.3870 

1.3663 

-4.0770 

-.2207 

.7500 

.7110 


1.0811 

-.6959 

2.8304 

-12.7121 

-.1615 

.8000 

.7532 

4.6583 

1.1049 

-1.0896 

6.4088 

-35.1357 

-.1015 

.8500 

.7945 

12.5662 

1.1517 

-2.0186 

17.4338 

-149-6938 

-.0535 

.9000 

.8349 

46.6378 

1.2275 

-4.3787 

65.3783 

-997.3960 

.0000 

1.0000 

.9129 

OO 

OO 

— OO 

OQ 

— OO 


k = -0.5 


I 


F(2 >) 

<J (1) (24) 

0 

0 

0.0007 

-0.0299 

0.0048 

-0.0602 

0.0105 

-0.0909 

0.0138 

-0.1115 

0.0187 

-0.1267 

0.1135 

-0.3009 

0.6430 

-0.703? 

3.6330 

-1.5871 

6.9763 

-2.2682 


q (4) 

E<°> 


r< 2 > 


.0000 


.0000 

.0000 

.0000 

.0002 

9.9600 

.0010 

-.0020 

.0000 

.0008 

4.9197 

.0030 

-.0103 

.0148 

.0064 

3.2141 

.0090 

-.0411 

.0670 


2.3371 

•0235 

-.1203 

.2307 


1.7947 

.0520 

-.3152 

.7600 

.5211 

1.5901 

.0755 

—5082 

1.4060 

1.3527 

1.4150 

.1088 

-.8277 

2.6949 

3.7802 

1.2620 

.1566 

-1.3796 

5.4353 

11.7558 

1.1254 

.2275 

-2.3946 

11.8168 

51.4676 





194.4562 

.8847 

.5248 

-9.0769 

81.6373 

954-643 

.7717 

.8743 

-21.9464 

306.6745 

16412.652 

.6555 

1.6804 

-71.6427 

1772.3621 


0000 


— OO 
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